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This  report  is  in  three  parts  corresponding  approximately  to  items 


of  the  contractual  work  statement. 
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P  Part  ljidescribes  calibration  of  the  radiometer  for  the  Background 
Measurements  and  Analysis  Program  (BMAP) .  ^ 
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Part  2  is  an  article  with  results  bearing  on  two-dim/lnsional  signal 


processing.  ^It  was  published  in 
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Part  3  consists  of  vugraphs  and  supplemental  text  for  two  BMAP 
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The  purpose  of 


1.1  PURPOSE  AND  SCOPE 


and  Analysis  Program  (BMAP7  is  to  col¬ 


lect  natural  infrared  background  data  for  use  in  simulating  the  performance  of  Infrared  Search 

and  Tracks  sets  (IRSTs).  Usually  the  objects  sought  by  an  1RST  are  effectively  point  sources; 

that  is,  their  images  are  unresolved  (no  larger  than  a  detector  element).  Thus,  the  point 

response  of  the  BMAP  sensor  and  electronics  must  be  well  characterized  for  the  intended 

data  use.  The  response  to  a  bar  source  must  also  be  determined  in  order  to  correct  the  data 

for -^certain  instrumental  artifact; described  infection  1.2.  An3~TiHIttj^  the  BMAP  instru- 
laJtsJD  * 

ments  should^be  radiometrically  calibrated  so  that  simulated  signals  added  to  the  background 

data  can  be  related,  at  least  approximately,  to  radiances  of  targets,  as  discussed  briefly  in 

references  1  and  2. 


The  point-  and  bar-source  calibrations  guarantee  that  the  data  truly  represent  the  spatial 

properties  of  the  background.  These  calibrations  are,  therefore,  more  important  than  the 

radiometric  calibration.  Radiometrically  uncalibrated  data  of  good  spatial  fidelity  can  be  used 

« 

very  effectively  for  comparative  studies  of  signal  processing  options,  as  references  1  and  2 
demonstrate.  Contrariwise,  radiometrically  calibrated  data  of  poor  spatial  fidelity  may  lead  to 
wrong  conclusions  if  used  for  simulations.  Nonetheless,  this  report  is  concerned  mainly  with 
radiometric  calibration  because,  overall,  it  is  the  most  difficult  of  the  three  calibrations  to 
implement,  perform  accurately  and  interpret.  Both  the  point-  and  bar-source  calibrations  are 
easy  to  implement  and  perform  accurately.  The  point-source  results  hardly  require  interpreta¬ 
tion,  although  the  apparatus  may  have  to  be  adjusted  to  obtain  the  proper  response.  The 
bar-source  results  are  hard  to  interpret  and  use  for  correcting  the  data,  but  R.  A.  Steinberg 
has  solved  the  problem  and  reduced  the  procedure  to  routine.3  For  these  reasons,  the  point- 
and  bar-source  calibrations  are  mentioned  in  this  report  only  to  give  a  complete  picture  of  the 
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calibration  procedure.  The  discussion  of  calibration  is  preceded  by  a  description  of  the 
apparatus,  which  gives  information  needed  for  understanding  the  calibration  and  using  the 
background  data. 

1.2  APPARATUS 

The  functional  parts  of  the  apparatus  are  two  optical  sensors,  their  analog  electronics, 

i 

and  a  data  acquisition  system.  This  report  concerned  mainly  with  the  optical  sensors  and 
analog  electronics,  which  make  up  the  radiometer.  Properties  of  the  apparatus  that  affect  cali¬ 
bration  and  use  of  the  background  data  are  given  in  Table  1  and  Figures  1-9. 

Figures  1  and  2  display  the  main  parts  of  the  optical  sensors  and  convey  an  idea  of  the 
radiometer’s  size.  The  analog  electronics  are  in  the  housing  with  the  optical  sensors.  In 
essence,  both  telescopes  are  permanently  focused  on  infinity;  the  longwave  (LW)  focus  can 
be  adjusted  slightly  to  correct  thermal  defocusing,  but  the  midwave  (MW)  optics  are  ther¬ 
mally  compensated  and  have  no  adjustment.  The  photosensitive  widths  and  heights  of  the 
detector  arrays  are  shown  larger  than  their  actual  scaled  sizes  in  Figures  1  and  2  in  order  to 
make  them  visible.  Figure  3  gives  the  sizes  and  layouts  of  the  detectors  to  scale.  The  tele¬ 
scopes’  blur  circles  (point-source  images)  are  about  (1/3)  miliiradian  in  diameter,  nearly 
matched  to  the  widths  of  the  detector  elements.  The  scan  mirror  oscillates  about  its  vertical 
axis  sweeping  images  of  the  scene  across  the  detectors.  Typical  point-source  outputs  are 
shown  in  Figure  4.  The  scan  period  is  0.5  seconds,  and  data  are  collected  during  a  single 
subinterval  of  each  period,  so  that  the  time  between  data  frames  is  0.5  seconds. 

The  scan  mirror  is  driven  by  a  triangular  wave.  The  periods  of  the  wave  and  mirror  are 
identical;  mirror  displacement  is  directly  proportional  to  wave  strength.  Consequently, 
increasing  the  wave  amplitude  increases  mirror  amplitude  and  scan  rate.  At  present  the  mir¬ 
ror  turns  approximately  ±3  (13/16)  degrees  about  its  vertical  axis  sweeping  the  line  of  sight 
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left  and  right  through  a  total  azimuthal  field  of  view  of  7(5/8)  degrees,  twice  the  mirror 
amplitude.  The  spatial  extent  of  the  data  interval  is  2.8125  degrees,  determined  by  a  shaft 
encoder;  the  temporal  extent  is  100  milliseconds,  determined  independently  as  an  electronic 
subdivision  of  the  driving  wave  period.  The  scan  rate  must,  therefore,  be  >28.125°/s  to 
keep  the  spatial  extent  of  the  data  interval  within  its  temporal  bounds. 

Figure  5  represents  the  radiometer’s  internal  sources  of  reference  radiation  for  one 
detector  array;  the  arrangement  for  the  other  is  the  same.  Cold  and  hot  references  are  pro¬ 
vided  by  motor-driven  mirrors  rotated  through  the  field  of  view.  The  cold  reference  mirror 
reflects  the  detector’s  own  radiation  back  to  it.  The  hot  reference  mirror  reflects  radiation 
from  an  internal  blackbody  onto  the  detector.  The  cold  reference  mirror  completely  fills  the 
field  of  view  for  about  2  milliseconds  to  allow  time  for  electronic  restoration  of  the  output 
level.  The  hot  reference  mirror  is  shaped  so  that  it  completely  fills  the  field  of  view  only  long 
enough  for  the  detector  to  respond  fully  to  the  internal  blackbody’s  radiation.  The  cold  refer¬ 
ence  mirror  passes  very  close  to  the  detector,  the  hot  reference  mirror  passes  farther  from  it. 
With  the  cold  feference  mirror  in  place,  the  detector  receives  hardly  any  radiation  except  that 
emitted  by  itself  and  the  mirror.  With  the  hot  reference  mirror  in  place,  the  detector  receives 
radiation  from  the  internal  blackbody  and  the  radiometer’s  interior.  When  the  detector  views 
the  scan  mirror  through  the  telescope,  it  receives  radiation  from  the  scene  and  stray  radiation 
from  outside  and  inside  the  radiometer. 

Figure  6  presents  block  diagrams  of  the  analog  electronics.  Each  detector  of  an  array  has 
its  own  electronics  like  those  shown.  The  amplifiers  and  cold  clamps  are  the  same  for  the 
two  arrays,  but  the  longwave  detectors  must  be  biased,  requiring  ac  coupling  to  eliminate  the 
bias  voltages  from  the  outputs.  This  potentially  complicates  the  radiometric  calibration,  but 
since  the  time  constants  of  the  ac  active  filters  are  about  30  times  the  data  interval,  the  effect 


of  ac  coupling  can  be  ignored  with  good  accuracy.  The  analog  outputs  go  to  the  data  acquisi¬ 
tion  system  where  they  are  sampled,  digitized,  and  recorded  on  magnetic  tape. 

The  output  format  is  displayed  in  Figure  7.  At  the  start  of  the  data  interval,  the  cold 
reference  mirror  rotates  into  place,  and  the  electronic  cold  clamp  restores  the  output  level  to 
ground  within  less  than  2.S%  of  the  full  scale  value,  10  volts  or  4096  arbitrary  units.  The 
output  then  rises,  peaks  and  falls  as  the  hot  reference  mirror  passes  through  and  out  of  the 
Held  of  view.  Afterward  the  detector  has  an  unobstructed  view  of  the  scene  for  the 
remainder  of  the  data  interval. 

Sampling  of  the  analog  outputs  is  controlled  by  the  scan  mirror  through  its  shaft  encoder 
which  takes  216  samples  in  360°,  i.e.,  512  samples  in  each  2.8125°  data  interval.  The  last  two 
samples  may  be  in  error,  however,  so  each  data  frame  contains  510  usable  values.  The  cold 
reference  occupies  the  first  10  values;  the  hot  reference  mirror  is  completely  out  of  the  field 
of  view  at  the  140th  sample;  the  remaining  370  samples  are  usable  scene  data  covering  2.03° 
or  35.5  milliradians  of  azimuth. 

The  scan  mirror  and  the  reference  mirrors  are  driven  by  waves  derived  from  the  same 
electronic  source,  an  IRIG  time  code  generator,  but  the  phases  of  the  reference-mirror  drives 
are  not  locked  to  the  scan-mirror  motion.  Consequently,  the  reference-mirror  rotation  can 
slip  relative  to  the  data  interval,  with  resultant  loss  of  the  cold,  or  cold  and  hot  references 
from  the  digital  output.  Failure  of  the  data  acquisition  system  can  also  result  in  loss  of  the 
references.  The  analog  outputs  always  contain  the  references,  but  not  always  during  the  time 
of  sampling  and  digitization. 

The  need  for  bar-source  calibration  results  from  pecularities  of  the  longwave  detector- 
array  and  its  electrical  connections  illustrated  in  Figure  8.  The  detectors  in  the  two  columns 
make  contact  with  a  common  gold  film  deposited  on  the  array  chip.  The  ends  of  the  gold 


film  are  connected  to  circuit  ground  by  two  platinum  wires  with  resistances  of  about  0.5  ohms 
each.  The  instantaneous  summed  current  I/„(f)  from  the  detectors  returns  to  the  gold  com¬ 
mon  through  the  platinum  wires.  The  resulting  voltage  across  the  wires  subtracts  from  the 
output  E„  ( t )  of  each  detector,  i.e., 

/„(/)  —  C„E„0)  —  D„T.In(t).  (1) 

Figure  9  illustrates  the  effect  of  this  resistive  coupling  on  the  output  produced  by  a  wide  bar 
source.  R.  A.  Steinberg  has  shown  how  to  determine  Cn  and  D„  from  measurements  of  x 
and  y  in  the  figure.3  Note  that  equation  (1)  is  valid,  and  the  coupling  constants  are  the  same, 
for  all  waveforms.  Use  of  the  bar-source  waveform  for  determining  C„  and  Dn  is  a  matter  of 
convenience.  That  is,  the  true  waveform  is  known  in  this  case  so  that  the  artifacts  are  easy  to 
identify.  In  addition,  the  equations  relating  the  coupling  constants  to  x  and  y  are  rather  sim¬ 
ple,  although  they  must  be  solved  by  numerical  iteration.  Once  Cn  and  D„  have  been  found 
from  the  bar-source  calibration,  however,  they  can  be  used  in  equation  (1)  with  measured 
outputs  /„(t)  from  any  source  to  find  the  true  outputs  E„(t )  from  that  source.  All  longwave 
data  are  corrected  in  this  way. 

1.3  RESULTS.  DISCUSSION,  CONCLUSIONS.  RECOMMENDATIONS 
1.3.1  Attempted  Calibration  at  Montauk  Point 

The  initial  attempt  at  radiometric  calibration  was  made  by  the  Naval  Surface  Weapons 
Center  (NSWC),  White  Oak,  during  background  measurements  at  Montauk  Point,  Long 
Island,  in  August  1983.  The  NSWC  procedure  was  to  scan  a  458C,  one  foot  square  blackbody 
placed  70  feet  from  the  radiometer. 

This  procedure  involves  a  number  of  difficulties  as  follows.  Corrections  must  be  made 
for  atmospheric  attenuation  of  the  blackbody  radiation.  Since  only  one  blackbody  tempera¬ 
ture  was  used,  the  offsets  (reference  radiation)  of  the  detectors  cannot  be  determined;  only 
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the  gains  (responsivities)  are  approximately  determinable.  Since  the  radiometer’s  telescopes 
are  focused  on  infinity,  the  blackbody  images  on  the  detectors  were  defocused  (Fig.  10).  The 
defocused  images  had  cores  of  unknown,  uniform  radiance  surrounded  by  borders  of  un¬ 
known,  nonuniform  radiance  (Fig.  11).  Since  the  telescopes  lie  above  and  below  the 
radiometer’s  centerline,  the  images  were  displaced  by  parallax  from  the  centers  of  the  detec¬ 
tor  arrays  (Figs.  11  and  12).  Consequently,  some  detectors  of  each  array  were  illuminated 
solely  by  the  nonuniform  borders  of  the  images  (Fig.  11).  (In  Fig.  11  the  detectors  are  num¬ 
bered  as  they  have  been  since  Montauk  Point.  Prior  to  and  at  Montauk  Point  they  were 
numbered  oppositely  with  1  at  low  and  16  at  high  elevation.  The  figures  of  Part  3 A  utilize 
the  old  numbering  system.) 


In  principle  the  radiances  of  the  defocused  images  could  be  calculated  and  used  to  derive 
the  detector  gains.  In  practice  this  wouid  be  very  difficult  to  do  and  would  give  results  of 
highly  questionable  accuracy.  It  is  recommended  that  the  Montauk  blackbody  scans  not  be 
further  analyzed  to  estimate  detector  gains  because  the  potential  returns  are  negligible  com¬ 
pared  to  the  effort  involved.  The  radiometer’s  spatial  response  was  well  calibrated  for  the 
Montauk  measurements,  so  that  the  background  data  are  valuable  without  radiometric  cali¬ 
bration,  as  explained  in  section  1.1. 

Several  possibilities  were  considered  for  changing  the  NSWC  procedure  to  overcome  the 
problems  of  defocusing  and  parallax.  The  obvious  solution  is  to  place  the  blackbody  farther 
from  the  radiometer.  However,  blackbody  size  and  distance  must  be  increased  proportionally 
to  keep  the  entire  detector  array  illuminated.  The  practical  limit  of  blackbody  size  makes  it 
impossible  to  increase  the  distance  sufficiently.  Defocusing  can  be  reduced  by  stops  (the 
pinhole-camera  effect)  or  eliminated  by  additional  lenses  ("spectacles"),  but  then  corrections 
must  be  made  for  the  effects  of  these  devices  in  order  to  obtain  a  calibration  for  the  radiome- 
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ter.  Obviously  little  is  accomplished  in  this  way.  In  short,  no  satisfactory  means  of  improving 
the  NSWC  procedure  has  been  found.  It  is  recommended,  therefore,  that  this  procedure  or 
variations  of  it  not  be  used  for  radiometric  calibration. 

1.3.2  Survey  of  Radiometric  Calibration 

Experience  with  the  NSWC  calibration  procedure  lead  to  a  review  of  the  literature  for 
the  purpose  of  developing  a  different  procedure.  References  4-8  were  the  most  useful 
sources  consulted.  Answers  to  three  questions  were  sought: 

What  is  a  radiometric  calibration? 

How  is  it  obtained? 

What  does  it  give  when  used  to  reduce  background  data? 

The  answers  to  these  questions  are  not  unique.  They  depend  on  the  data  usage,  the  instru¬ 
ment,  the  scene,  and  the  field  conditions. 

As  already  noted,  BMAP  data  are  to  be  used  for  comparing  and  evaluating  IRST  designs 
by  simulation  of  their  performance.  This  use  requires  data  in  radiance  units  for  simulating 
optical  sensors  with  arbitrary  apertures  and  instantaneous  fields  of  view  (IFOV).  The  BMAP 
radiometer  should,  therefore,  be  calibrated  directly  in  radiance  units.  Calibrations  in  units  of 
irradiance,  intensity,  or  power  would  have  to  be  converted  to  radiance  units.  This  would 
require  measurements  of  the  radiometer’s  aperture  and/or  IFOV,  which  would  involve  addi¬ 
tional  effort  and  error. 

A  radiometer  responds  to  radiation  in  a  broad  passband  of  wavelengths  from  Xj  to  X2. 
The  output  depends  on  received  power,  the  product  of  received  radiance  with  the  constant 
IFOV  and  aperture  area.  Thus,  output  can  be  expressed  as  a  function  of  spectral  radiance 
L(\)  integrated  over  the  passband.  For  data  reduction,  however,  the  inverse  relation  giving 
integrated  radiance  as  a  function  of  output  v  is  wanted,  i.e., 


R0  fk*  r(\)  L(\)  d\- f(v),  (2) 

where  Rq  r(\)  is  the  radiometer’s  spectra!  response.  R0  is  the  instrument’s  responsivity  at 
the  arbitrary  wavelength  \0  where  the  response  function  r(\)  is  normalized  to  unity. 

Radiometric  calibration  is  the  process  of  determining  integrated  radiance  as  a  function  of 
radiometer  output  from  measurements  of  the  outputs  produced  by  a  standard  source.  The 
radiance  of  the  standard  source  should  be  varied  over  a  range  such  that  the  outputs  encom¬ 
pass  the  range  of  outputs  produced  by  the  unknown  sources.  Calibrations  always  depend  on 
the  radiometer’s  spectral  response;  measured  radiances  of  unknown  sources  do  also  if  equa¬ 
tion  (2)  is  used  for  calibration.  However,  the  desired  quantities  are  the  integrated  radiances 
of  the  unknowns  without  instrumental  dependence.  There  are  two  ways  of  calibrating  so  that 
the  measured  radiances  are  independent  of  the  radiometer’s  spectral  response.  Neither 
involves  equation  (2),  and  both  require  matching  the  spectrum  of  the  standard  to  that  of  the 
unknown  at  some  point  along  the  line  of  sight. 

The  best  of  the  two  special  calibrations  gives  a  relation  between  radiometer  output  and 
integrated  clutter  radiance,  i.e., 

£  Le  d\  «■  /j(v),  (3) 

where  Lc  is  the  clutter’s  spectral  radiance.  To  determine  this  relation  experimentally,  it  is 
necessary  to  use  a  standard  source  with  the  same  spectrum  as  the  clutter  and  to  place  the 
source  beside  the  clutter  so  that  the  atmosphere  affects  source  and  clutter  radiation  identi¬ 
cally.  BMAP  conditions  forbid  doing  either  because  the  scenes  contain  objects,  including 
clouds,  at  different  ranges  with  different  spectra.  Thus,  the  desired  calibration  is  unattainable 
in  BMAP  because  (1)  the  necessary  standard  sources  are  unavailable,  and  (2)  they  could  not 
be  deployed  properly  is  they  were  available.  (This  type  of  calibration  is  best  because  the 
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measured  clutter  radiances  do  not  depend  on  either  the  radiometer’s  spectral  response  or 
meteorologic  conditions.  The  calibration  depends  on  both  but  includes  them  in  such  a  way 
that  its  use  to  reduce  the  background  data  makes  the  results  independent  of  these  conditions. 
It  is  not  even  necessary  that  the  instrument's  passband  include  the  integration  range  wholely 
or  partly;  it  is  only  necessary  that  radiation  in  the  two  wavebands  be  uniquely  related.) 

The  second  special  calibration  gives  a  relation  between  radiometer  output  and  integrated 
scene  radiance  at  the  radiometer’s  entrance  aperture,  i.e., 

[t a(d)Le  +  La(d))d\  -  /2(v),  (4) 

where  d  is  range  to  the  clutter  and  (rfl  La)  are  the  spectral  (transmittance,  radiance)  of  the 
intervening  atmosphere.  A  standard  source  with  the  same  spectrum  as  the  scene  radiation  at 
the  entrance  aperture  is  needed  for  this  calibration.  Such  a  source  is  unavailable  for  BMAP 
because  the  scene  spectra  at  the  entrance  aperture  are  unknown  and  variable.  (This  type  of 
calibration  gives  measured  radiances  independent  of  the  radiometer’s  spectral  response  but 
not  of  meteorologic  conditions.  Matching  the  spectra  of  the  standard  and  unknown  sources 
always  provides  independence  from  spectral  response  but  not  from  meteorologic  conditions 
when  the  spectrum  includes  atmospheric  effects.) 

Since  the  spectra  of  the  standard  source  and  scene  cannot  be  matched  under  BMAP  con¬ 
ditions,  it  is  necessary  to  use  an  arbitrary  standard.  The  usual  choice  is  a  blackbody.  The 
measured  scene  radiances  then  depend  on  both  the  radiometer’s  spectral  response  and 
meteorologic  conditions,  which  must  be  recorded  when  the  background  data  are  collected. 
The  BMAP  radiometer’s  view  port  is  a  convenient  fiducial  point  for  radiance  regardless  of  the 
blackbody’s  location.  The  calibration  equation  in  this  case  is 

R°l  r[TaLb  +  La]d\  -  /3(v), 


(5) 


bration  should  not  be  used  outside  its  range.  The  same  types  of  results  for  all  midwave  and 
longwave  channels  are  given  in  Table  5  and  6,  where  the  variances  are  the  squared  standard 
errors  of  estimate,  and  the  95%  prediction  limits  are  calculated  at  values  of  u  for  which  the 
output  v  is  about  25,  50,  75,  and  100%  of  full  scale. 

« 

Further  measures  of  calibration  quality  can  be  obtained  as  follows.  Assume  ideal  condi¬ 
tions,  that  is,  errorless  calibrations  and  measurements  made  with  an  extended  source  per¬ 
fectly  uniform  in  space  and  time.  Then  radiances  calculated  from  the  dc  coupled  midwave 
outputs  should  be  identical  across  and  along  the  channels  at  each  and  every  sample  number. 
Radiances  calculated  from  the  ac  coupled  longwave  outputs  should  be  identical  across  the 
channels  at  each  sample  number  but  should  droop  slightly  with  increasing  sample  number. 
This  can  be  visualized  with  the  aid  of  Figure  15  by  imagining  the  quantity  on  the  vertical  axes 
to  be  radiance.  Beyond  the  140th  sample  and  under  ideal  conditions,  the  midwave  plots 
should  be  in  a  horizontal  plane;  the  longwave  plots  should  be  in  a  plane  sloping  downward  to 
the  right,  ignoring  slight  curvature  of  the  droop.  Furthermore,  under  ideal  conditions, 
midwave  (or  longwave)  radiances  at  the  same  sample  number  should  be  the  same  in  all 
frames.  Departures  from  ideal  conditions  cause  differences  in  radiances  at  the  same  sample 
number  both  within  and  between  frames. 

W,  t  actually  is  observed  in  contrast  to  the  ideal?  Figures  20  and  21  show  cross¬ 
channel-average  radiance  calculated  at  each  sample  number  from  data  recorded  with  a  near 
extended  blackbody  at  50.0°C.  Figures  22  and  23  show  the  cross-channel  standard  deviations 
of  the  radiances  at  each  sample  number.  Under  ideal  conditions,  the  midwave  averages 
should  be  on  a  horizontal  line;  the  longwave  averages  should  be  on  a  line  sloping  down  to  the 
light,  again  ignoring  slight  curvature  of  the  droop;  and  the  standard  deviations  should  be 
zero.  The  departures  from  ideal  behavior  are  small  and  apparently  random  except  for  a  slight 
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1.3.8  Radiometric  Calibrations  Derived  and  Evaluated 

The  radiometric  calibration  on  tape  BF7  was  reduced  as  follows.  First,  the  radiometer 
outputs  were  processed  to  remove  "transients,"  features  which  rise  or  fall  more  rapidly  than 
the  radiometer’s  point  response  allows.  Next,  resistive  coupling  artifacts  (section  1.2)  were 
removed  from  the  longwave  outputs.  (These  corrections  are  applied  routinely  to  all  data 
R.  A.  Steinberg  originated  both.)  Then  values  of  the  response-weighted  radiance  integral 
<rLb>  were  calculated  for  each  external  blackbody  temperature  in  Table  3  using  the 
response  functions  of  Figures  16  and  17.  Finally,  equation  (26)  was  fit  by  least  squares  to  the 
fractional  outputs  u  and  integrated  radiances  <rLb>  to  determine  offsets  A  and  gains  B  for 
all  midwave  and  longwave  channels.  (Longwave  data  at  a  blackbody  temperature  of  30°C 
were  not  used  because  the  outputs  had  no  cold  references.) 

Results  for  one  midwave  and  one  longwave  channel  are  presented  in  Figures  18  and  19. 
The  equations  fit  the  data  so  well  that  the  plots  fail  to  show  differences  between  the  two. 
These  are  revealed  only  by  the  statistical  measures  of  the  fits.  The  standard  error  of  estimate 
is  the  standard  deviation  between  the  observed  radiances  and  those  calculated  from  the  cali¬ 
bration  equations.  The  sample  correlation  is  the  fraction  of  the  relation  between  u  and 
<rLb>  accounted  for  by  the  linear  dependence.  With  reference  to  Figure  18,  for  example, 
the  95%  confidence  intervals  and  limits  of  prediction  have  the  following  meanings:  The  pro¬ 
bability  is  95%  that  the  offset  A  lies  in  the  range,  —0.25 1  ±3.4  x  10-3  W/m2-sr.  If  the  output 
is  one-quarter  of  full  scale  (f.s.),  the  probability  is  95%  that  the  true  integrated  radiance  lies 
within  ±6.0  x  10~3  W/m2sr  of  the  value  given  by  the  calibration  equation.  Of  course  these 
measures  refer  only  to  the  precision  of  this  one  calibration,  not  to  its  reproducibility  or  accu¬ 
racy.  Notice  that  the  limits  of  prediction  are  smallest  at  the  center  of  the  calibration  range 
and  increase  with  the  magnitude  of  the  departure  from  the  center;  this  is  one  reason  a  cali- 
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gain  or  to  subtract  a  constant  voltage  from  the  analog  output.  Reducing  gain  is  inadvisable 
since  the  longwave  clutter  lies  mainly  in  the  top  half  of  the  present  digital  range,  so  that  actu¬ 
ally  an  increase  of  gain  is  desirable.  Subtracting  a  strictly  constant  voltage  from  the  analog 
output  is  awkward  and  difficult.  The  best  solution  is  to  use  an  A/D  converter  with  a  max¬ 
imum  greater  than  the  present  4096,  which  also  would  allow  an  increase  of  gain. 

1.3.7  Field  Measurements 

During  the  week  of  10-14  September  1984,  cloud  background  data  were  collected  from  a 
rooftop  site  at  Raytheon  Missile  Systems  Division  overlooking  Hanscom  Air  Force  Base. 
After  equipment  setup  on  the  first  day,  data  were  collected  in  good  weather  on  the  second  to 
fourth  days,  but  not  on  the  fifth  day  when  the  weather  was  bad.  A  radiometric  calibration 
was  put  at  the  end  of  each  data  tape  according  to  the  procedure  for  field  calibration  given  in 
section  1.3.5. 

Eight  tapes  were  made  in  this  way.  The  tapes  and  supporting  documentation  (notes, 
radiometer  logs,  meteorologic  logs,  and  camera  logs)  are  already  in  the  sponsor’s  possession. 
Plots  of  selected  data  frames  were  made  to  spot  check  the  tapes.  The  plots  show  massive 
failures  of  the  data  acquisition  system  with  the  result  that  only  the  seventh  and  eighth'tapes 
are  useful.  Tables  3  and  4  summarize  the  radiometer  logs  for  these  tapes,  showing  the  time 
required  for  the  calibrations  and  how  they  are  put  on  the  tapes.  Either  the  midwave  or  the 
longwave  calibration  takes  less  than  five  minutes.  The  two  had  to  be  made  separately  because 
the  blackbody  was  not  large  enough  to  cover  both  telescopes.  With  a  commercially  available 
blackbody  twice  as  high,  the  two  calibrations  could  be  made  simultaneously  in  less  than  five 
minutes  by  recording  the  outputs  in  the  A  (alternating)  mode.  From  the  standpoint  of  the 
time  required,  the  procedure  employed  is  certainly  practicable.  It  is  even  feasible  to  calibrate 
more  often  than  once  per  tape  should  the  calibrations  be  temporally  or  thermally  erratic. 


30°C.  However,  the  resistive  coupling  artifacts  in  Figure  14  were  not  completely  removed  by 
coupling  parameters  derived  about  a  year  earlier  from  outputs  of  a  wider  single-bar  source 
(Figure  9).  Clearly  the  coupling  parameters  changed  with  time.  Accordingly,  new  coupling 
parameters  were  determined  from  the  more  recent  measurements,  even  though  x  and  y  of 
Figure  9  are  not  as  easily  measured  from  a  narrow  output.  For  this  reason  a  wider  bar  target 
should  be  used  in  the  future. 

Radiometric  calibrations  made  in  the  laboratory  have  not  yet  been  reduced,  although  the 
detector  outputs  have  been  plotted  and  assessed.  (The  196  plots  have  been  given  to  the 
sponsor  separately  to  avoid  making  this  report  too  bulky.  Figure  15  displays  examples  from 
the  set.)  The  plots  show  that  the  cold  reference  level  is  sometimes  missing,  owing  either  to 
failure  of  the  data  acquisition  system  or  to  slippage  between  the  data  interval  and  the  refer¬ 
ence  mirror  rotation.  For  this  reason,  the  entire  midwave  calibration  at  a  radiometer  tem¬ 
perature  of  30°C,  and  certain  frames  of  the  calibrations  at  25  and  35°C,  may  be  unusable.  In 
addition,  stray  radiation  tends  to  distort  or  saturate  the  hot  reference  peak  at  high  radiometer 

and  standard  source  temperatures.  For  this  reason,  the  midwave  calibration  with  the  radiom- 

* 

eter  at  40°C  and  the  standard  source  at  50°C  may  be  unusable.  Despite  these  shortcomings 
prospects  are  good  for  successfully  reducing  most  of  the  midwave  calibrations.  The  opposite 
is  true  of  the  longwave  calibrations.  Those  at  radiometer  temperatures  of  25  and  30°C  were 
made  but  not  recorded  because  of  operator  error.  At  35  and  40°C  the  cold  reference  is  miss¬ 
ing,  and  the  hot  reference  is  often  saturated  by  stray  radiation.  Of  all  the  longwave  laboratory 
calibrations,  only  the  one  at  a  radiometer  temperature  of  20°C  appears  usable. 

Saturation  of  the  hot  reference  by  stray  radiation  probably  makes  the  longwave  radiome¬ 
ter  unusable  above  30°C  —  86°F.  This  serious  limitation  should  be  removed.  Saturation 
occurs  in  the  digital,  not  the  analog  output.  Potential  solutions  are  to  reduce  the  amplifier 


3.  Near  extended  source  scans 

a.  Place  source  over  view  port  in  contact  with  radiometer  housing. 

b.  Record  four  frames  of  data  from  each  detector  array. 

c.  Do  b  with  source  at  10,  20,  30,  40,  50°C. 

4.  Do  1-4  with  radiometer  at  20,  25,  30,  35,  40°C. 

The  following  simpler  procedure  is  used  in  the  field: 

1.  Point  source  scans 

a.  Center  image  from  portable  collimated  source  vertically  on  pixel. 

b.  Record  four  frames  of  data. 

c.  Do  a-b  for  channels  2,  8,  15  of  detector  array. 

d.  Do  a-c  for  both  arrays. 

2.  No  bar  source  scans. 

3.  Near  extended  source  scans 

a.  Samf  ptocedure  as  in  laboratory. 

b.  Record  scans  at  end  of  each  data  tape. 

1.3.6  Laboratory  Measurements 

The  BMAP  radiometer  was  calibrated  in  the  laboratory  at  Raytheon  Missile  Systems 
Division,  Bedford,  Massachusetts  on  23  August  1984  by  the  procedure  in  the  previous  sec¬ 
tion.  The  notes,  logs,  and  digital  tapes  made  at  that  time  have  been  turned  over  to  the  spon¬ 
sor. 

Typical  point-source  outputs  with  the  radiometer  at  20.6°C  are  shown  in  Figure  4.  Fig¬ 
ure  14  is  a  representative  bar-source  output  at  the  same  radiometer  temperature.  The  shapes 
and  widths  of  these  outputs  indicate  an  acceptable  spatial  response  for  the  radiometer  at  20- 


temperature.  For  that  reason  the  responses  must  be  calibrated  in  both  the  laboratory  and  the 
field.  Laboratory  calibrations  are  needed  to  provide  baseline  measures  of  the  responses.  In 
the  laboratory,  conditions  can  be  controlled  and  reproduced;  measurements  can  be  made 
slowly  with  care  and  repeated  to  be  sure  they  are  correct.  Then  a  simpler,  faster  calibration  in 
the  field  is  sufficient  to  verify  the  radiometer’s  response  by  comparison  with  the  laboratory 
results. 

The  laboratory  calibration  setup  is  illustrated  in  Figure  13.  Point  and  bar  sources  are 
placed  at  infinity,  in  effect,  through  use  of  a  collimating  mirror.  The  radiometer  is  put  on  a 
micrometer  stage,  allowing  the  point-source  image  to  be  centered  carefully  on  each  detector 
for  determining  the  point  responses.  The  bar  source  is  high  enough  to  illuminate  all  detec¬ 
tors  of  an  array  simultaneously.  For  radiometric  calibration  an  extended  biackbody  source  is 
placed  over  the  radiometer’s  view  port  in  contact  with  the  instrument’s  housing.  The 
radiometer’s  temperature  is  varied  systematically  to  obtain  the  temperature  dependence  of  the 
point,  bar  and  radiometric  responses.  To  investigate  the  temporal  stability  of  the  responses, 
all  of  these  measurements  are  made  at  intervals  of  days,  weeks  or  months. 

The  step-by-step  procedure  for  laboratory  calibration  is  as  follows: 

1.  Point  source  scans 

a.  Center  point-source  image  vertically  on  pixel. 

b.  Record  four  frames  of  data. 

c.  Do  a-b  for  all  channels  of  detector  array. 

d.  Do  a-c  for  both  arrays. 

2.  Bar  source  scans 

a.  Center  5 0°C  bar  target  vertically  on  detector  array. 

b.  Record  two  frames  of  data  from  each  array. 
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The  radiometer  output  with  the  standard  source  in  place  (ignoring  atmospheric  radiance)  is 

V  -  y</'(rLb  +  /)>  +  (y  </>  -1-  8).  (23) 

Subtracting  equation  (20)  from  equation  (23),  dividing  the  result  by  equation  (22),  and 


rearranging  terms  gives 


<f'rLb>  + 


</'(/-  /o)> 


</'('*-/<>»  '  *  </'(4 -4>»  ’ 


where  as  before  u  —  (v  —  v0)/(vA  -  v0)  is  the  fractional  output  with  respect  to  the  hot  and 
cold  reference  voltages.  Solving  equation  (24)  for  <f'rLb>  gives  the  broadband  calibration 
equation  for  the  BMAP  radiometer, 

< /'t Lb >  -  <f'(lh  -  /(,)>  u  -  </'(/  -  /0)>.  (25) 

For  convenience,  equation  (25)  is  written  in  terms  of  a  response  function  r  and  broad¬ 
band  offset  and  gain  parameters  A  and  B. 

<rLb>  —  A  +  Bu  (26) 

The  response  function  includes  the  detector’s  response  /'  and  all  radiative  losses  r  in  the 
atmosphere  and  optics  between  source  and  detector.  Emission  from  the  optics  is  included  in 
the  offset  through  the  stray  radiance  /.  The  offset  and  gain  include  the  response  of  the  detec¬ 
tor  but  not  the  electronics.  With  broadband  radiation,  expressing  the  output  by  the  ratio  u 
compensates  the  calibration  parameters  for  electronics  drift  but  not  detector  drift.  Nonethe¬ 
less,  radiometer  temperature  probably  will  affect  the  broadband  parameters  as  it  does  the 
monochromatic  ones,  i.e.,  gain  probably  will  increase  slowly  and  offset  more  rapidly  with  tem¬ 
perature. 


1.3.5  Calibration  Procedure 


The  spatial  and  radiometric  responses  of  the  BMAP  radiometer  may  vary  with  time  and 


most  of  lh,  so  the  increase  of  (lh  —  /0)  with  temperature  should  be  slow.  This  will  result  in  a 
slow  increase  of  the  gain  (lk  —  l0)/r  with  temperature.  The  difference  (/  —  /„)  will  increase 
with  temperature  more  rapidly  than  (/*  —  /0),  since  /  consists  entirely  of  stray  radiation  from 
the  radiometer’s  main  optics  and  interior.  The  offset  (/  -  /q)/ t  should,  therefore,  increase 
with  temperature  more  rapidly  than  the  gain. 

B.  Broadband  Radiation 

Integration  over  wavelength  changes  these  results  very  little,  as  will  now  be  demon¬ 
strated.  For  convenience,  the  integration  will  be  performed  over  all  wavelengths  so  that  the 
limits  are  effectively  set  by  the  cut-offs  of  the  spectra  and  responses.  This  avoids  conceptual 
problems  arising  from  the  fact  that  scene  radiation,  standard  radiation,  internal  reference  radi¬ 
ation  and  stray  radiation  follow  different  optical  paths  and  have  different  spectra  at  the  detec¬ 
tor.  The  results  are  the  same  as  those  obtained  by  integrating  over  the  spectral  passband  and 
accounting  separately  for  the  differences  in  paths  and  spectra,  provided  the  passband  is 
chosen  as  described  in  section  1.3.2.  For  simplicity,  the  integration  will  be  denoted  by  angular 
brackets,  i.e., 

<*>  -  fo  * 

From  equation  (12)  the  spectrally  integrated  radiometer  output  is 

v -></'/>+ (y</>  +  8).  (19) 

With  the  cold  or  hot  reference  mirror  in  place,  the  outputs  are 

v0  -  7 <A>>  +  (y</>  +  8),  (20) 

v*  -  y<A>  +  (y</>  +8).  (21) 

The  amplification  factor,  which  now  does  not  include  the  detector  response  /',  is  given  by 


v*  -  v0-  </'(/*  -  /0)  >  y. 


(22) 


for  stray  radiation  and  r  for  losses  of  radiation  between  the  source  and  detector.  (Atmos¬ 
pheric  radiance  is  ignored.)  With  the  standard  source  in  view,  the  voltage-radiance  relation  is 


v  —  y/'(rLb  +  /)  +  (yf  +  8) 

-  yf'rLb  +  ( yf'l  +  yf  +  8). 


(16a) 

(16b) 


Subtracting  equation  (13)  from  equation  (16b)  and  dividing  the  result  by  equation  (15) 


V/i  ~  v0  I  ib  i0 


‘h  *0 


This  relation  expresses  the  radiometer’s  output  as  a  fraction  u  of  the  output  voltage  with  the 
hot  reference  mirror  in  place,  all  voltages  being  measured  from  the  cold  reference  level  v0. 
Solving  equation  (17)  for  Lb  gives  a  monochromatic  calibration  equation  for  the  BMAP 


radiometer. 


l„-lo 


Equation  (18)  does  not  involve  the  responses  of  the  detectors  and  electronics.  These  are 
accounted  for  in  the  variable  u  through  v0  and  vA,  which  are  measured  before  each  scan. 
Thus,  the  calibration  is  compensated  for  drifts  of  the  detectors  and  electronics  over  times 
longer  than  the  scan  time. 


The  bracketed  terms  in  equation  (18)  correspond  to  the  radiometer  gain  and  offset 
parameters.  This  allows  the  temperature  dependence  of  the  parameters  to  be  predicted.  The 
loss  factor  t  probably  depends  weakly  on  radiometer  temperature.  The  radiance  /0  increases 
slowly  as  the  temperature  of  the  cold  reference  mirror  increases.  The  difference  (lh  —  l0)  will 
increase  with  temperature  since  lh  includes  stray  radiation  from  the  internal  reference  optics 
and  the  radiometer’s  interior.  However,  radiation  from  the  internal  hot  source  accounts  for 


and  add  a  voltage  8  (Fig.  6),  so  that  the  voltage-radiance  relation  for  the  whole  radiometer 
becomes 

v  -  yf'l  +  (yf  +  8).  .  (12) 

(Throughout  this  section,  the  same  symbol  v  is  used  for  all  voltage  outputs  with  the  precise 
meaning  indicated  by  context  and  subscripts.)  The  factors  y  and  6  are  constant  only  if  the 
electrical  and  thermal  variables  of  the  system  are  fixed.  The  added  voltage  S  is  chosen  so  that 
the  radiometer  output  is  v0  when  the  detector  is  completely  covered  by  the  cold  reference 
mirror  (Figs.  5  and  7).  During  the  time  of  complete  coverage,  the  radiance  /0  at  the  detector 
comes  from  emission  by  the  mirror  and  reflection  of  the  detector’s  own  radiation  back  to  it. 
Hardly  any  stray  radiation  reaches  the  detector  during  this  time.  The  voltage-radiance  relation 
with  the  cold  reference  mirror  in  place  is 

v0  —  yf'lo  +  (y/  +  8).  (13) 

The  hot  reference  mirror  follows  the  cold  reference  mirror  through  the  field  of  view. 
With  the  hot  reference  mirror  in  place,  most  of  the  radiance  lh  at  the  detector  comes  from 
the  internal  hot  source'and  the  optics  which  transfer  its  radiation  to  the  detector.  However, 
stray  radiation  from  the  radiometer’s  interior  also  reaches  to  the  detector  at  this  time  (Fig.  5). 
The  peak  hot  reference  voltage  (Fig.  7)  is 

v„-rA  +  (y/  +  5).  (14) 

The  difference  between  the  hot  and  cold  reference  voltages  gives  a  relation  for  the  amplifica¬ 
tion  factor  y /'. 

v*  ~  v0  -  (/A  -  l0)yf  (15) 

For  calibration  the  radiometer’s  view  port  is  covered  by  an  extended  blackbody  of  radi¬ 
ance  Lb.  Under  these  conditions  the  radiance  at  the  detector  is  ( rLb  +  /),  where  /  accounts 


actual  use  an  equation  which  accounts  for  specific  features  of  the  BMAP  radiometer  is 
needed.  That  equation  is  obtained  by  analyzing  the  response  of  the  detectors  and  electronics 
to  scene,  internal  reference,  and  stray  radiation.  The  analysis  also  indicates  how  the  calibra¬ 
tion  parameters  should  depend  on  the  radiometer  temperature.  Only  the  midwave 
radiometer’s  response  is  analyzed  since,  as  already  explained,  the  longwave  radiometer 
responds  in  approximately  the  same  way  because  of  its  very  long  electronic  time  constant. 
The  response  to  monochromatic  radiation  will  be  derived  first,  then  the  analysis  will  be 
extended  to  broadband  radiation.  This  will  avoid  obscuring  the  essential  features  of  the 
response  in  the  fog  of  integration  over  the  spectral  passband. 

A.  Monochromatic  Response 

Assume  that  the  detector  output  voltage  v  is  an  analytic  function  of  the  optical  (O), 
thermal  (77,  and  electrical  (£)  variables  of  the  instrument  and  environment, 
v  -  HiO,T,E).  If  all  variables  except  power  p  received  by  the  detector  are  fixed,  v  —  hip), 
and  since  power  is  the  product  of  throughput  (IFOV  x  aperture  area)  with  radiance  /  at  the 
detector,  v  —  /(/).  This  relation  can  be  expanded  in  a  Maclaurin  series. 

v  -  /( 0)  +  /'(0)/  +  y/"(0)/2  + - .  (10) 

The  coefficients  of  the  series  are  functions  of  all  the  system  variables  except  l  and,  there¬ 
fore,  are  constants  only  if  all  variables  except  radiance  at  the  detector  are  fixed.  It  will  be 
assumed  that  the  first  two  terms  of  the  series  adequately  represent  the  voltage-radiance  rela¬ 
tion.  Then 

V-/  +  /7,  (11) 

where  the  point  of  series  expansion  has  been  omitted  for  simplicity. 

The  midwave  detector  output  is  dc  coupled  to  electronics  which  amplify  it  by  a  factor  y 


In  summary,  this  section  has  (1)  demonstrated  that  the  BMAP  radiometer  should  be 
calibrated  with  a  blackbody  radiance  standard,  (2)  described  the  rational  basis  of  such  a  cali¬ 
bration,  and  (3)  explained  the  meaning  of  scene  radiances  determined  with  the  calibration. 

1.3.3  Selection  of  Blackbody  Calibration  Technique 

The  possible  ways  of  calibrating  a  radiometer  with  a  blackbody  standard  can  be  classified 
very  simply:  The  source  can  be  large  or  small,  and  it  can  be  placed  close  to  the  radiometer  or 
far  from  it  (Table  2).  The  most  important  attributes  of  a  calibration  technique  for  BMAP 
usage  are  listed  at  the  left  of  Table  2.  Some  of  these  derive  from  the  fact  that  the  technique 
must  be  usable  in  the  field  as  well  as  the  laboratory.  A  few  comments  will  clarify  the  mean¬ 
ing  of  and  need  for  the  listed  attributes.  A  technique  which  does  not  require  focusing  on  the 
standard  is  best  since  the  radiometer’s  telescopes  provide  little  or  no  focusing.  Corrections 
for  atmospheric  attenuation  of  the  standard  radiation  can  be  avoided  by  placing  the  source 
close  to  the  radiometer.  Background  radiation  should  not  reach  the  detector  during  calibra¬ 
tion.  Radiation  from  the  standard  should  pass  through  all  parts  of  the  aperture  as  the  back¬ 
ground  radiation  does  when  data  are  collected.  The  standard  should  be  commercially  avail¬ 
able  to  avoid  the  time  and  expense  of  developing  a  special  source.  For  use  in  the  field,  the 
calibration  procedure  must  be  fast  and  simple.  The  source  must  be  easy  to  handle,  control 
and  position  in  the  field  of  view,  and  it  must  illuminate  the  whole  array  to  allow  simultaneous 
recording  of  all  detector  outputs.  Table  2  indicates  advantage,  neutrality,  disadvantage,  or 
uncertainty  for  the  possible  calibration  techniques  with  respect  to  the  listed  attributes. 
Clearly,  the  near-extended-source  technique  is  best  for  calibration  of  the  BMAP  radiometer. 

1.3.4  BMAP  Calibration  Equation 

Equation  (6)  of  section  1.3.2  exhibits  the  calibration  equation  in  general  form,  but  for 
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where  the  parameters  g  and  g'  are  determined  from  measurements  of  v  and 

r[raLb  +  La]dk  -  JQ  r[raLb  +  La]dk.  (9) 

Now  the  limits  of  integration  in  equation  (8)  include  all  values  of  the  integrand,  and  so  equa¬ 
tion  (8)  describes  the  radiometer’s  response  to  either  standard  or  scene  radiation  with 
L  —  [raLb  +  La)  or  L  —  Wa(d)Lc  +  La(d)].  That  is  not  necessarily  true  of  equation  (6) 
because  the  calibration  passband  may  exclude  significant  values  of  r[Ta(d)Lc  +  La(d )]  if  the 
scene  radiation  is  very  strong  outside  the  passband.  In  other  words,  if  two  sources  produce 
the  same  radiometer  output  v, 

X  rL,d\  -  fQ  rL2dk 

always,  but 

rL\d\  *  fK*rL2d\ 

unless  the  limits  of  integration  include  all  significant  parts  of  both  integrands.  Generally  the 

passband  cannot  be  chosen  to  satisfy  this  condition  because  the  scene  spectra  are  unknown 

and  variable.  The  alternative  is  to  choose  the  passband  for  calibration  to  satisfy  equation  (9). 

« 

In  that  case  scene  radiances  calculated  from  equation  (8)  are  values  of 

r[ra{d)Lc  +  La(d)]d\, 

that  is,  response-weighted  scene  radiance  at  the  view  port  integrated  over  all  wavelengths. 

With  the  spectral  response  written  R0r(\)>  the  response  shape  r(X)  must  be  unity  at 
wavelengths  where  the  radiometer’s  responsivity  is  Rq.  The  spectral  response  is  normalized  . 
by  selecting  a  wavelength  X0  at  which  r(X0)  is  made  to  equal  one.  The  choice  of  X<j  is  arbi¬ 
trary  and  determines  values  of  both  R o  and  r(X).  Customarily  Xo  is  chosen  to  be  the 
wavelength  at  the  peak  of  the  spectral  response.  Both  the  spectral  response  and  the  measured 
radiances  are  then  described  as  "peak  normalized"  or  "nomalized  to  unity  at  the  peak." 
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where  Lb  is  blackbody  spectral  radiance  and  (rc,  La)  are  the  spectral  (transmittance,  radi¬ 
ance)  of  the  atmosphere  between  the  blackbody  and  the  view  port.  The  bracketed  term  is 
spectral  radiance  at  the  view  port  during  calibration,  and  the  response  Function  includes  all 
radiative  losses  between  the  view  port  and  the  detector.  (Any  position  between  the  view  port 
and  the  blackbody  could  serve  as  the  fiducial  point  with  appropriate  changes  in  the  response 
and  radiance  terms.)  Ordinarily  equation  (5)  is  divided  by  R0  and  only  two  or  three  terms  in 
the  series  expansion  of  g(v)  =■  f3  (v)/R0  are  retained,  with  the  result 

rLd\  =-  g(v) 

—  #(0)  +  #'(0)v  +  -j£"(0)v2  H - .  (6) 

The  spectral  radiance  at  the  view  port  is  denoted  by  L  in  equation  (6)  so  that  it  can  represent 
either  standard  or  scene  radiation  in  the  ensuring  discussion  of  passband  determination  and 
the  meaning  of  measured  scene  radiance. 

In  practice  the  radiometer  responds,  however  weakly,  to  radiation  of  all  wavelengths. 
Limiting  the  integration  in  equation  (6)  to  a  finite  passband  is  an  approximation.  The  exact 
form  of  equation  (6)  is 

f  1  rLd\  +  f  1  rLd\  +  f  rLd\  —  #(v).  (7) 

*/0  Jk2 

The  first  and  third  integrals  must  be  negligible  compared  to  the  second  for  equation  (6)  to 
approximate  the  exact  equation  closely.  This  requires  the  passband  to  include  all  significant 
values  and  all  long  weak  "tails"  of  the  integrand  rL.  The  limits  K\  and  X2  are  chosen  to  satisfy 
this  requirement  during  calibration;  hence,  the  interval  Xt  to  X2  should  be  called  the  "calibra¬ 
tion  passband."  When  the  calibration  satifies  the  stated  requirement,  equation  (6)  can  be 
written 


Jo  rLd\  -  g(0)  +  #'(0)v, 


(8) 
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curvature  of  the  longwave  averages  in  the  wrong  direction.  The  greatest  departure  of  the 
midwave  averages  from  a  horizontal  line  through  them  is  about  0.06%.  The  greatest  depar¬ 
ture  of  the  longwave  averages  from  a  smooth  curve  through  them  is  about  0.03%,  and  the 
total  droop  is  about  0.1%.  The  midwave  and  longwave  standard  deviations  are  about  0.1% 
and  0.2%  of  the  respective  average  radiances.  These  results  indicate  high  uniformity  of  the 
standard  source  over  the  scanned  part  of  its  surface.  They  also  demonstrate  excellent  internal 
consistency  of  the  calibrations;  there  are  no  large  systematic  errors  between  channels.  The 
smallness  of  the  droop  justifies  use  of  the  midwave  calibration  equation  (26)  for  the  longwave 
calibration  as  well. 

1.3.9  Interpretation  of  Measured  Radiances 

The  passband  limits,  Aj  and  X2,  are  chosen  so  that  blackbody  radiation  outside  the 
passband  contributes  insignificantly  to  the  integral  of  response-weighted  blackbody  spectral 
radiance  Lb,  i.e., 

J~  rLbd\  -  rL„dX.  (27) 

The  equation  for  the  radiometer’s  response  as  determined  by  calibration  using  equation  (26) 
is  then 

J  rLd\  —  A  +  Bu,  (28) 

where  L  is  either  blackbody  or  scene  spectral  radiance  at  the  instrument’s  view  port.  The 
radiometer’s  response  function  r  includes  all  radiative  losses  from  the  view  port  to  the  detec¬ 
tor  and  is  normalized  to  unity  at  the  response  peak  (Figs.  16  and  17).  With  a  cluttered  scene 
in  view 

L  —  Ta(d)Lc  +  La(d),  (29) 

where  Lc  is  clutter  spectral  radiance,  d  is  range  to  the  clutter,  and  (rfl,La)  are  the  spectral 


(transmittance,  radiance)  of  the  intervening  atmosphere.  Thus  the  measured  scene  radiances 
are  values  of 

J^~  r[ra(d)Lc  +  La(d)]d\, 

that  is,  response-weighted  scene  radiance  at  the  view  port  integrated  over  all  wavelengths. 

The  measured  quantities  are  not  those  desired,  which  are  the  spectral  or  integrated 
clutter  radiances.  These  can  be  derived  from  the  measured  quantities  as  follows.  Write  the 
scene  spectral  radiance  in  the  form, 

rfl(X,  d)  LC(X)  +  La(k,  d)  «•*  Los(x),  (30) 

where  L0  is  a  constant  and  s(X)  is  the  scene  spectrum  at  the  view  port.  If  both  r(X)  and 
s(X)  are  known,  L0  and  the  scene  spectral  radiance  at  the  view  port  can  be  determined  since 

L0  rsdk  -  j£“  r[ra{d)Lc  +  La(d)]d\.  (31) 

The  definite  integrals  are  simply  numbers,  those  on  the  right  hand  side  being  the  measured 
scene  radiances.  Now  if  range  to  the  clutter  and  meteorologic  conditions  also  are  known, 
clutter  spectral  radiance-  can  be  calculated  from  equation  (30)  in  the  form, 

Z,e(X)  -  (I0s(X)  -  La  (X,  d)]/ra  (X,  d),  (32) 

using  values  of  atmospheric  radiance  and  transmittance  given  by  a  model  such  as  Lowtran  6. 

Presently  scene  radiances  at  the  view  port  cannot  be  determined  from  the  response- 
weighted  BMAP  radiances  because  scene  spectra  at  the  view  port  are  not  known.  Clutter 
radiances  cannot  be  determined  for  the  same  reason  and  also  because  ranges  to  the  clutter  are 
not  known.  In  addition,  the  response  functions  of  Figures  16  and  17  may  be  inadequate  for 
determining  scene  spectral  radiances;  if  the  scene  spectra  are  very  strong  outside  the 
passbands,  values  of  the  responses  which  are  insignificant  in  equation  (27)  may  be  needed  to 
solve  equation  (31)  for  £0.  This  seems  improbable  under  BMAP  conditions,  but  the 


possiblity  should  be  kept  in  mind.  The  quantities  displayed  in  Figures  16  and  17  are  called 
passbands  and  response  functions  for  calibration,  because  they  were  determined  specifically 
for  calibrating  the  radiometer  and  may  be  unsuitable  for  other  uses.  It  is  recommended  that 
scene  spectra  and  ranges  to  the  clutter  be  measured  if  possible,  and  that  the  response  func¬ 
tions  be  determined  more  accurately  outside  the  calibration  passbands  if  necessary. 
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1.5  FIGURES 


1.  Radiometer  Top  View  (1/2  Scale):  Scan  Mirror,  Telescope,  Detector 

2.  Radiometer  Side  View  (1/2  Scale):  Scan  Mirror,  Telescopes, 
Detectors 

3.  Detector-Array  Geometry 

A.  Typical  Point-Source  Outputs 

5.  Radiometer  Internal  References 

6.  Midwave  and  Longwave  Analog  Electronics  Diagrams 

7.  Output  Format 

8.  Origins  of  Longwave  Resistive-Coupling  Effect 

9.  Longwave  Bar-Source  Outputs  with  Resistive-Coupling  Artifacts 

10.  Blackbody  Image  Formation  during  Attempted  Radiometric  Calibration 
at  Montauk  Point 

11.  Inverted,  Defocused,  Parallactic,  Blackbody  Image  Crossing  Midwave 
Detector  at  Montauk  Point 

12.  Parallactic  Effects  in  Montauk  Calibration  Setup 

13.  Top  View  of  Laboratory  Calibration  Setup 

1A,  Longwave  Output  from  a  Multiple  Narrow  Bar  Source 

15.  Typical  Frames  of  Radiometric  Calibration  Data 

16.  Midwave  Response  Function  for  Calibration 

17.  Longwave  Response  Function  for  Calibration 

18.  Midwave  Calibration,  Channel  6,  Tape  BF7 

19.  Longwave  Calibration,  Channel  6,  Tape  BF7 

20.  Cross-Chainnel  Averages  of  Midwave  Radiances 

21.  Cross-Channel  Averages  of  Longwave  Radiances 

22.  Cross-Channel  Standard  Deviations  of  Midwave  Radiances 

23.  Cross-Channel  Standard  Deviations  of  Longwave  Radiances 
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OMETER 


RADIOMETER  TOP  VIEW  (1/2  SCALE):  SCAN  MIRROR,  TELESCOPE,  DETECTOR 
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Internal  hot  reference  voltage 
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FIGURE 


Mr)  (Arbitrary  Units)  *2(t)  (Arbitrary  Units) 


sample  number 


LONGWAVE  BAR-SOURCE  OUTPUTS  WITH  RESISTIVE-COUPLING  ARTIFACTS 


l.UUU  X  1U  “  W/M 
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CHANNEL  AVERAGES  OF  MIDWAVE  RADIANCES 


FIGURE  17 


FIGURE  16 


RRD IOMETER 
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T  -  20. 6C,  BRR  T  -  50. OC,  RRW  DRTR 

DAYS  12  HRS  6  MIN  54.8  SEC 


FRAME  :  166 

CHANNEL  s  A 
8AND  :  LONGWAVE 


5.  81ackbodj  sou-ree 

3.  EOI  modal  H5  153, 

4,7  mil  aperture. 

b.  EOI  model  P  K06  Si 

with  6  bar  target 

c.  EOI  model  Pl406Sl 

w/o  i  bar  target 


1.  Optical  bench 

l.  Movable  stage, 

horizontal  4  vertical 
micrometer  adjustments 

3.  Radiometer 

4.  Off-axiS  paraboloid, 

10 14  dia.j  100"  i.  1. 


TOP  VIEW  OF  LABORATORY  CALIBRATION  SETUP 


FIGURE  13 


Detector 


PARALLACTIC  EFFECTS  IN  MONTAUK  CALIBRATION  SETUP 
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INVERTED,  DEFOCUSED,  PARALLACTIC,  BLACKBODY  IMAGE 
CROSSING  MIDWAVE  DETECTOR  AT  MONTAUK  POINT 


FIGURE  11 
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BLACKBODY  IMAGE  FORMATION 

ATTEMPTED  RADIOMETRIC  CALIBRATION  AT  MONTAUK  POINT 


1.6  TABLES 

1.  System  Parameters  for  Data  Recording 

2.  Blackbody  Calibration  Options  and  BMAP  Selection  Criteria 

3.  Summary  of  Radiometer  Log  for  Tape  BF7 

4.  Summary  of  Radiometer  Log  for  Tape  BF8 

5.  Midwave  Calibration  Parameters  and  Statistics,  Tape  BF7 

6.  Longwave  Calibration  Parameters  and  Statistics,  Tape  BF7 
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TABLE  1:  SYSTEM  PARAMETERS  FOR  LATA  RECORDING 


Item 


Midwave  Spectral  Band 
( 3 . 8-5 . 0  urn) 


Longwave  Spectral  Band 
(7.3-11.9  pm) 


Telescope 

4  inch  diameter  aperture 

6  inch  focal  length 

Same  as  midwave 

Detector  array 

InSb  photovoltaic,  77°E 

16  detectors  in  2  columns  of  8a 
See  Figure  1  for  layout 

HgCdTe  photoconductive,  77°K 
Same  as  midwave*3 

II 

Field  of  view  (mR) 

instantaneous 

total 

0.37  az.  x  0.40  el.c 

49.09  az.  x  5.40  el.d 

0.33  az.  x  0.33  el. 

49.09  az.  X  5.33  el. 

Spectral  response 

See  Figure  16 

See  Figure  17 

3dB  electronic  passband  (Hz) 

0  -  1000 

0.3  -  1000e 

Point-source  output 

full  width  at  half  maximum 

See  Figure  2  for  shape 

0.37  mRf 

Same  as  midwave 

0.33  mR 

Recorded  random  noise 

2.5  mv 

Same  as. midwave 

NEI  (I/cm2) 

2  x  10'14 

2  x  10"13 

extent 

7(5/8)°  of  scene 

Same  as  midwave 

data  interval 

2.8125°,  92.2  ms 

•i 

period 

0.5  s  (time  between  data  ints.) 

it 

rate  during  data  Interval 

30.5°/s  *  532  mR/sK 

•i 

dwell  time 

0.695  msb 

0.620  msh 

Sampling  time 

interval 
intervals /dwell 
data  interval 
cold  reference 
cold  &  hot  references 
usable  scene  data 


<<  dwell  time 
0. 09606  mR 
3. 851 

512  samples 
10  samples^ 
140  samples 
370  samples 


Same  as  midwave 


3.44 

Same  as  midwave 


Digitization 


0-10  v  "  12  bits  *  0-4096  decimal  Same  as  midwave 


(a)  Detectors  14  and  16  dead  at  Montauk  Point 

(b)  Detector  1  dead 

Values  at  Montauk  Point:  (c)  0.33  az.  x  0.33  el.,  (d)  49.09  az.  x  5.33  el.,  (e)  0.5  -  1000, 

(f)  0.33  mR,  (g)  36.0°/s  “  628  mR/s,  (b)  0.525  ms,  (i)  3.44,  (J)  23 


4-  :  Advantage  0  :  Neutral  -  :  Disadvantage 

Soarce  Location  Distant 

and  Size 

Delect  ion 

Criteria  Small  Large  Small  Large 


Soarce  Location 
and  Size 


Focus  not  required 

— 

— 

+■ 

4- 

Atmospheric  attenuation 

— 

— 

4- 

4- 

Stray  background  radiation 

— 

4- 

—  j 

4- 

Aperture  filled 

+ 

4- 

— 

4- 

Direct  radiance  responsivitg 

- 

4- 

4- 

4* 

Source  availability 

+■ 

2 

Z 

-h 

Whole  array  illuminated 

— 

2 

+ 

+  ’ 

Ease  of  handling  4-  Control 

- 

— 

0 

4* 

Positioning  not  critical 

0 

0 

““ 

4- 

BLACKBODY  CALIBRATION  OPTIONS  ANI 

TABLE  2 

1  BMAP  SEI 

_ECTI0N  CRITERIA 

PLACE:  Bedford,  Massachusetts,  USA  DIGITAL  TAPE:  BF7 

42°28,42‘*Ni  71°17'24',W,  250  ft.  above  mean  sea  level 
DATE:  13  September  1984  VIDEO  TAPE:  45664302 

SUMMARY  OF  RADIOMETER  LOG 

Object  Scanned  Local  Time  Record  Frame  Line  of  Sight*  Tempera ture(°C) 


Ho. 

Description 

of  Start 

Mode* 

Nos . 

Azimuth 

Elevat ' 

'n  Optics 

IntBB 

ExtBB 

SCE33 

Hole  In  cloud  deck 

08:52:21.3 

A  ' 

001-061 

108 

12 

25.2 

37.4 

— 

SCE34 

•• 

08 : 55 : 36  .*8 

" 

062-163 

112 

15 

25.1 

It 

— 

SCE3S 

•• 

08:58:39.8 

ft 

164-265 

109 

19 

•• 

tt 

— 

SCE36 

09:07:51.3 

« 

266-347 

130 

11 

24.9 

•• 

— 

SCE37 

M 

09:11:47.8 

8/8 

348-389 

113 

09 

24.2 

tf 

— 

SCE38 

•• 

09:18:07.8 

A 

390-461 

093 

08 

23.3 

tt 

— 

XBB9 

External  black  body 

09:24:42.3 

MW 

462-465 

— 

— 

23.0 

37.4 

50.0 

»» 

It 

09:25:34.8 

•• 

466-469 

— 

— 

22.9 

tf 

40.0 

M 

ft 

09:26:30.8 

" 

470-473 

— 

— 

ft 

M 

30.0 

fl 

»• 

09:27:28.8 

•« 

474-477 

— 

— 

ft 

ft 

20.0 

M 

•t 

09:28:25.8 

tt 

478-481 

— 

H 

- 

10.0 

XBB9 

External  black  body 

09:29:31.8 

LI 

482-485 

— 

— 

22.9 

37.4 

10.0 

M 

»« 

09:30:37.8 

" 

486-489 

— 

— 

22.7 

tt 

20.0 

M 

•• 

09:32:07.3 

tf 

490-493 

. - 

— 

22.6 

tt 

30.0 

n 

»# 

09:33:08.3 

ft 

494-497 

— 

— 

If 

tt 

40.0 

tf 

ft 

09:34:05.3 

t« 

498-501 

— 

— 

#• 

tt 

50.0 

CC7 

Closed  cover 

09:34:54.3 

LI 

502-505 

— 

— 

22.6 

37.4 

— 

tt 

ft 

09:35:51.3 

HI 

506-509 

— 

— 

t* 

•' 

— 

SCE39 

Hole  In  cloud  deck 

09:38:33.3 

8/8 

510-551 

120 

14 

22.6 

37 

— 

*A  »  alternate  midwave  and  longwave  frames.  MW  or  LW  **  midwave  or  longwave  output  only. 

8/8  ■  longwave  outputs  (5-12)  in  data  channels  (1-8)  plus  midwave  outputs  (5-12)  in  data 
channels  (9-16), 

Azimuth  Id  degrees  clockwise  from  geographic  N.  Elevation  in  degrees  aboye  horizontal. 
IntBB ,  ExtBB  -  Internal,  External  Black  Body 


SUMMARY  OF  RADIOMETER  LOG  FOR  TAPE  BF7 


TABLE  3 


56 


PLACE:  Bedford,  Massachusetts.  USA  DIGITAL  TAPE:  BF8 

42°28,42"M,  71°17‘24"W,  250  ft.  above  mean  sea  level 
DATE:  13  September  1984  VIDEO  TAPE:  45664302 


SUMMARY  OF  RADIOMETER  LOG 


No. 

Object  Scanned 

Description 

Local  Time 
of  Start 

Record 

Mode* 

Frame 
Nos . 

Line  of 
Azimuth 

Sight* 
Elevat 'n 

Temperature(°C) 
Optics  IntBB  ExtBB 

SCE40 

Stratocumulus 

15:36:55.3 

A 

001-059 

233 

33 

29.6 

37.7 

— 

SCE41 

.15:42:34.3 

8/8 

060-101 

242 

30 

27.2 

37.6 

— -- 

SCE42 

" 

15:46:15.8 

•• 

102-143 

253 

19 

27.3 

ft 

— 

SCE43 

Thin,  rapidly  moving 
clouds 

15:48:16.8 

" 

144-265 

229 

22 

27.4 

M 

— 

SCE44 

Stratocumulus 

15:58:11.8 

266-307 

275 

24 

27.8 

** 

— 

SCE45 

» 

16:19:31.3 

308-369 

238 

06 

26.9 

re 

— 

SCE46 

•• 

16:21:31.3 

370-411 

248 

08 

26.6 

rr 

— 

XBB10 

External  black  body 

16:31:58.8 

MV 

412-415 

— 

— 

26.4 

37.6 

50.0 

t* 

t« 

16:33:11.8 

•• 

416-419 

— 

— 

•• 

<» 

40.0 

•  • 

** 

16:34:30.3 

420-423 

— 

26.5 

ee 

30.0 

M 

‘ 

16:35:24.3 

424-427 

— 

— 

N 

•« 

20.0 

f* 

i* 

16:36:41.3 

" 

428-431 

— 

— 

26.4 

•e 

10.0 

XBBIO 

External  black  body 

16:37:27.8 

LW 

432-435 

— 

— 

" 

10.0 

•  • 

»• 

16:38:43.3 

»» 

436-439 

— 

— 

26.3 

re 

20.0 

»» 

»» 

16:39:54.8 

•» 

440-443 

— 

— 

•« 

•r 

30.0 

ft 

M 

16:41:03.8 

444-447 

— 

— 

26.2 

" 

40.0 

" 

»• 

16:42:05.8 

" 

448-451 

— 

— 

•> 

»r 

50.0 

CCS 

Closed  cover 

16:44:50.3 

LW 

452-455 

— 

— 

26.1 

37.6 

_ 

M 

16:45:21.8 

MV 

456-459 

— 

— 

n 

» 

_ 

SCE47 

Stratocumulus 

16:59:03.3 

A 

460-550 

257 

06 

25.2 

H 

_ 

*A  “  alternate  midwave  and  longwave  frames.  MW  or  LW  —  midwave  or  longwave  output  only. 

®/8  "  longwave  outputs  (5-12)  in  data  channels  (1-8)  plus  midwave  outputs  (5-12)  in  data 
channels  (9-16). 

Azimuth  In  degrees  clockwise  from  geographic  N.  Elevation  in  degrees  above  horizontal. 
IntBB.  ExtBB  ”  Internal,  External  Black  Body 


SUMMARY  OF  RADIOMETER  LOG  FOR  TAPE  BF8 


TABLE  A 


I 

I 

CHANNEL  1 

PARAMETERS 

A  1  B 

1  1 
I  VAR IANCE  | 

95%  CONFIDENCE 
A  1 

INTERVALS  1 

B  t 

1  1 

-0. 195 

1  2.  23 

1  3.  591E-06  1 

2,  319E-03  1 

2. 927E-03  1 

3  1 

-0.  201 

1  2.  27 

1  3.  812E-06  1 

2.  397E-03  1 

3. 060E-03  ! 

3  1 

-0. 232 

1  2.  30 

1  3.  194E-06  | 

2.  843E-03  1 

3. 620E-03  1 

4  ( 

-0. 231 

1  2.  32 

1  3.  720E-06  | 

2.  405E-03  1 

3. 093E-03  1 

S  1 

-O. 238 

t  2.  33 

1  3.  237E-06  1 

2. 275E-03  1 

2, 920E-03  1 

6  1 

-0. 231 

1  2.  36 

1  7.  109E-06  t 

3. 339E-03  1 

4. 361E-03  1 

7  1 

-0.  268 

1  2.  38 

1  1.  205E-06  1 

1.  394E-03  1 

1. B03E-03  1 

8  1 

-0.  236 

1  2.  38 

1  4.  916E-06  I 

2.  B00E-03  I 

3, 631E-03  1 

9  1 

-0.  276 

1  2.  41 

1  3  141E-07  | 

9. 148E-04  I 

1.  197E-03  1 

10  1 

-0. 267 

I  2.  43 

1  3.  771E-06  1 

3.  031E-03  t 

4. 037E-03  ! 

11  : 

-O.  273 

!  2.  42 

i  5.  943E-06  1 

3.  105E-03  1 

4,  08 IE -03  1 

13  1 

-0.  230 

1  2.  39 

i  4.  965E-06  1 

2.  B03E-03  1 

3. 6B8E-03  I 

13  1 

-0.  263 

S  2.  40 

1  4,  407E-06  | 

2  661E-03  I 

3. 4B7E-03  ! 

14  I 

-0.  244 

1  2.  39 

!  4.  624E-06  1 

2  699E-03  1 

3, 537E-03  1 

15  1 

-0. 231 

1  2.  39 

i  4.  132E-06  I 

2  567E-03  1 

3. 367E-03  ! 

16  1 

-O.  227 

1  2.  37 

!  3.  707E-06  1 

2  972E-03  1 

3. 920E-03  i 

95%  PREDICTION  LIMITS 


!  I 

:  CHANNEL  1 

0.  300  1 

U 

0.  630  : 

0  960  1 

1 

1  30  1 

/  1  1 

4.  271 £-03  1 

4,  091  £-03  1 

4  !3-'E-03  : 

4  406E-03  1 

i  2  : 

4.  399E-03  i 

4  213E-03  ! 

4  263E-03  1 

4  559E-03  1 

i  3  ! 

5.  145E-03  ! 

4.  919E-03  1 

4  97BS-03  ! 

3  326E-03  1 

:  4  i 

4.  35 IE— 03  1 

4  1 6 1 £-03  1 

4  218E-03  1 

4  524E-03  l 

:  5  : 

4.  066E-03  1 

3  883E-03  1 

3. 934E-03  1 

4  223E-03  1 

i  6  : 

6.  019E-03  ; 

3. 752E-03  1 

3, 839E-03  i 

6  282E-03  1 

!  7  i 

2.  480E— 03  I 

2. 368E-03  1 

2. 403E-03  ! 

2.  586E-03  1 

i  B  1 

3.  004E-O3  ! 

4  783E-03  1 

4. B39E-03  1 

3.  234E-03  1 

!  9  I 

1 . 620E-03  ! 

1. 347E-03  1 

1  572E-03  1 

1.  697E-03  1 

:  to  ! 

3.  42  1  £—03  1 

3.  180E-03  1 

3.  275E-03  1 

3  707E-03  1 

:  ui 

3.  503E-03  1 

3.  238E-03  1 

3  349E-03 

5  779E-03  1 

1  12  1 

3.  023E-03  1 

4.  B03E-03  ! 

4  8B8E-03  1 

3  274E-03  I 

1  13  1 

4.  739E-03  ! 

4.  32BE-03  1 

4.  604E-03  ! 

4  96BE-03  1 

1  14  ! 

4.  047E-O3  1 

4  637E-03  1 

4.  7  1 BE-03  1 

3  092E-03  1 

1  13  1 

4.  396E-03  1 

4.  395E-03  i 

4  469E-03  1 

4  B20E— 03  1 

:  16  i 

5.  379E-03  ! 

3  1 3 1 £-03  1 

5  24  IE-03  ! 

3  652E-03  1 

CORRELATION:  l.  00 


MIDWAVE  CALIBRATION  PARAMETERS  AND  STATISTICS,  TAPE  BF7 


TABLE  5 


1 

I 

i 

CHANNEL  ! 

PARAMETERS 

A  1  B 

• 

i  VARIANCE 

1 

1 

95%  CONFIDENCE  INTERVALS 

A  1  B 

4 

t 

I 

2  1 

-12.  9 

1  5b,  2 

1  5. 109E-03 

1 

0.  218 

1  0.  284 

• 

• 

1 

3  t 

-12.7 

1  56.  1 

!  5. 241E-03 

1 

0.  220 

1  0. 287 

• 

1 

4  t 

-12.  7 

1  55.  9 

t  4.  030E-03 

» 

i 

0.  193 

1  0.  251 

1 

» 

• 

5  ! 

-12.  8 

i  56.  2 

I  3.  B48E-03 

i 

0.  189 

t  0.  246 

• 

• 

1 

i 

6  1 

-12.  8 

1  55.  8 

t  3. 849E-03 

i 

0.  189 

t  0.  245 

« 

• 

i 

7  i 

-12.  5 

!  55.  6 

!  4. 435E-03 

» 

i 

0.  201 

1  0.  262 

1 

■ 

t 

« 

8  : 

-12.  6 

5  55.  4 

i  4. 327E-03 

i 

0.  199 

1  0.  257 

« 

• 

« 

i 

9  1 

-12.  B 

1  55.  5 

S  5.  699E-03 

[ 

0.  230 

1  0.  296 

1 

i 

io  : 

-13.  0 

(  55.  4 

!  4. 379E-03 

• 

• 

0.  202 

1  0.  259 

• 

• 

t 

■ 

11  1 

-12.  8 

i  55.  0 

i  4. 17BE-03 

J 

0.  197 

t  0.251 

1 

1 

i 

12  1 

-12.  B 

1  54.  5 

1  3.  277E-03 

i 

0.  174 

1  0.  220 

1 

i 

13  1 

-13.  3 

!  54.  8 

i  5. 668E-03 

• 

0.  232 

1  0.  291 

: 

i 

14  1 

-13.  2 

!  54.  0 

!  3. 852E-03 

i 

» 

0.  191 

1  0. 237 

i 

• 

i 

15  I 

-14.  1 

!  54.  8 

!  5. 705E-03 

i 

0.  237 

1  0.  292 

f 

i 

: 

16  1 

-13.  7 

!  53.  9 

I  4. 439E-03 

! 

0  207 

1  0.  254 

• 

95%  PREDICTION  LIMITS 


i 

i 

1 

CHANNEL  I 

0.  250 

1  0.  500 

U 

1 

0.  750 

1  1.  00 

I 

I 

i 

2  1 

0.  213 

1  0.  174 

1 

0.  158 

1  0.  173 

1 

i 

3  : 

0.  216 

1  0.  176 

1 

0.  160 

I  0.  175 

1 

i 

4  1 

0.  189 

!  0.  154 

1 

1 

0.  140 

I  0. 153 

1 

1 

5  1 

O.  185 

1  0.  151 

1 

0.  137 

1  0.  150 

1 

i 

6  1 

0.  185 

1  0.  151 

! 

0.  137 

1  0.  149 

1 

1 

i 

7  1 

0.  198 

i  0.  162 

f 

0.  147 

1  0.  160 

1 

i 

8  1 

0.  196 

!  0. 160 

i 

0.  145 

1  0.  158 

1 

» 

i 

9  1 

0.  226 

i  0.  184 

1 

t 

0.  167 

1  0.  181 

1 

» 

( 

10  1 

O.  199 

!  0.  162 

I 

0.  146 

1  0.  158 

1 

i 

11  1 

0.  194 

:  0.  158 

I 

0.  143 

1  0  154 

I 

! 

12  : 

0.  172 

1  0.  141 

1 

0.  127 

I  0.  136 

( 

1 

13  : 

O.  228 

:  0.  186 

1 

0.  167 

1  0.  178 

1 

1 

14  | 

0.  188 

i  0.  154 

1 

t 

0.  138 

I  0.  146 

1 

1 

« 

15  i 

0.  231 

i  0.  188 

• 

1 

0.  168 

i  0.  177 

1 

» 

16  1 

0.  204 

:  0  166 

♦ 

« 

0  148 

1  0.  155 

» 

CORRELATION:  1.00 


LONGWAVE  CALIBRATION  PARAMETERS  AND  STATISTICS,  TAPE  BF7 


TABLE  6 


SYNCHRONOUS  DEDICATED  IROCEDIRES 


Activation 


1.  Power  on  equipment. 

2.  When  connected,  identify  your  terminal  to  the  network  by  entering 

D$$SON  XXXXXX  'enter  your  sign-on  identification  code  for  XXXXXX) 
(Example:  >$$SON  AA1234)  &  K  0  1  /  £ 

3.  The  system  will  respond  with: 

UNIVAC's  Tel con  Network-Level  XXX  DCPID-XXX 
SESSION  PATH  OPEN 

4.  Identify  your  site  to  the  DROLS  system  by  transmitting  the  following: 


DSGNONS/TERMINAL  ID  _  On  *1  OaL  £  / A-  77  3^  </ 

(Example:  DSGNONS/ABCDE)  ^  '  '  ' 

5.  The  system  will  respond  with: 

*MSG  ONI  SIGN-ON  ACCEPTED 

6.  Enter  your  terminal  ID  as  in  the  past.  Entering  of  the  terminal  ID  at  this 
point  will  be  eliminated  in  the  near  future. 

7.  DROLS  commands  remain  the  same. 


Termination 

1.  Enter  @TniM@  and  transmit. 

2.  The  system  will  respond  with: 

THIS  TERMINAL  HAS  BEEN  TERMINATED 
*MSG  D07  PLEASE  SIGN  OFF  TERMINAL 

3.  Disconnect  (sign  off)  from  the  network  by  transmitting: 
C>$$SOFF 

4.  The  system  will  respond  with: 

SESSION  PATH  CLOSED 


^  . 


5.  Power  off  equipment. 


TERMINAL  USER  CONDITION  MESSAGE  DESCRIPTION 


TERMINAL  USER  CONDITION  MESSAGE  DESCRIPTION 


PART  2 


SIMULATION  OF  MID-INFRARED  CLUTTER  REJECTION. 
THRESHOLD-SENSOR-SIZE  EFFECTS  WITH  LMS-FILTERED  NOISE. 


Published  in  Applied  Optics,  15  February  1985 

by 

Longmire,  F.  D.  Bryant,  J.  D.  Wilkey,  and  A.  F.  Milton 


04  3  0  c?79 


Abstract 


The  relation  between  an  adaptive  threshold  sensor's  size  and  its 

performance  is  analyzed  numerically  using  LMS-filtered  noise.  The 

analysis  links  signal  detection  to  size  through  the  threshold 

parameter,  sampling  variance,  and  a  clutter-edge  effect.  A  principle 

is  established  for  determining  by  statistical  analysis  of  natural 

background  data  whether  a  2-D  threshold  sensor  may  be  advantageous. 

It  is  not  under  the  conditions  of  this  simulation.  Critical, 

optimum,  and  best  sizes  are  given  for  a  1-D  threshold  sensor 

following  a  1-D  LMS  filter  with  input  from  an  optical  sensor  having 

-13  2 

NEI  (noise  equivalent  irradiance)  *=  1.5  x  10  W/cm  and  spatial 
resolution  0.15  x  0.36  mR  (milliradians ) . 


-2- 


I.  Introduction 


Wide-field  infrared  (IR)  surveillance  sets  must  detect  signals 
against  backgrounds  with  uniformly  radiant  parts  (often  clear  sky) 
and  cluttered  parts  of  variable  radiance  (often  sunlit  clouds). 
Techniques  for  detecting  signals  in  uniform-sky  noise  are  well 
developed,  but  those  for  detecting  signals  in  clutter  noise  -  or 
simply  rejecting  clutter  to  avoifi  false  alarms  -  are  not.  This  state 
of  affairs  exists  chiefly  because  theoretical  analysis  using  a  simple 
noise  model  gives  the  optimum  method  of  detecting  signals  in  uniform 
sky,  but  not  in  clutter.  With  clutter  an  alternative  approach  is  to 
compare  signal  processing  options  using  a  computer,  recorded  natural 
background  noise,  and  programmed  algorithms  to  simulate  the 
processing.  This  approach  identifies  the  best  option  without 
requiring  each  to  be  tested  in  a  surveillance  set.  Simulation  does 
not  give  theoretically  optimum  solutions,  but  that  is  not  perforce  a 
disadvantage.  For  example,  the  theoretically  optimum  method  of 
detecting  point-source  signals  in  clear'  sky  fails  altogether  in 
clutter  and  has,  therefore,  limited  value  for  wide-field  IR 
surveillance.  In  this  case,  a  compromise  between  performance  in  clear 
sky  and  clutter  may  be  more  useful  than  optimum  performance  in 
either.  A  possible  compromise  was  studied  and  several  signal 
processing  options  were  compared  in  the  article  preceding  this  one. 

Both  articles  explore  clutter  rejection  and  signal  detection  by 
electronic  filters  and  'hreahold  sensors  operating  on  the  output  of  a 


mid-IR  (4. 0-4. 8  pm)  optical  sensor  scanned  over  clear  sky  and  clouds 
with  very  bright,  backlit  edges.  Detailed  descriptions  of  the  data, 
the  algorithms  and  the  method  of  simulation  appear  in  the  previous 
article.^  (An  appendix  abstracting  those  descriptions  is  provided 
here  at  the  request  of  the  editor- and  reviewers  to  make  this  article 
reasonably  self-contained.)  The  previous  article  compared  the 
performance  of  two  threshold  sensors  -  one  fixed,  the  other  adaptive 
-  following  several  one-dimensio.nal  ( 1  — D)  linear  filters.  The 
combination  of  a  1-D  adaptive  threshold  sensor  and  a  1-D  LMS  spatial 
filter  was  found  to  suppress  extreme  cloud  clutter  effectively  and 
to  detect  simulated  signals  almost  as  well  in  clutter  as  in  clear 
sky.  This  combination  merits  further  study,  since  the  performance  of 
an  adaptive  threshold  sensor  depends  on  its  size  -  its  spatial  extent 
or  the  number  of  data  used  to  estimate  the  local  noise  statistics  - 
but  only  one  size  was  studied  in  the  previous  article.  Here  we 
report  the  results  of  further  work  to  determine  how  changing  the 
size  of  a  1-D  threshold  sensor  following  a  1-D  LMS  filter  affects 
signal  detection. 

The  relation  between  an  adaptive  threshold  sensor's  size  and  its 
performance  has  already  been  investigated  theoretically,  e.g.,  see 
references  2-7.  Most  of  this  work  has  dealt  with  processing  of  radar 
or  sonar  signals,  but  the  principles  apply  also  to  processing  of 
infrared  signals.  The  theory  is  not  directly  applicable  to  real 
cases,  however,  because  it  must  employ  low  fidelity  clutter  models 
and  other  simplifying  assumptions  to  render  the  problem 
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false  alarm  permitted  in  all  the  data  always  occurs  in  the  clutter. 


Except  for  a  clutter-edge  effect  discussed  below,  the  threshold 
parameters  vary  with  frame  size  as  the  sampling  s.d.s  do.  This  means 
that  there  are  certain  critical  sample  sizes  beyond  which  the  thresh¬ 
old  parameter  and  signal  detection  improve  very  slowly.  The  critical 
sample  size  depends  on  the  fractional  rate  of  decrease  of  the 
sampling  s.d.,  which  is  obtained  by  differentiating  equation  (2). 


1 

d(  ) 

n 

(c r„/a) 

dN 

n 


(3) 


If  £  is  the  fractional  rate  of  decrease  at  the  critical  sample  size 

N  ,  then  N  ■  a/2£.  From  Figs.  4A  and  4B  the  critical  sample  size 
c  c 

in  uniform  sky  appears  to  be  about  100  uncorrelated  noise  values  (250 

data),  which  agrees  with  the  theoretical  result  for  a  radar  set 

operating  in  a  constant  noise  environment,  reference  2,  p.  421.  This 

-3  -2 

gives  £  «  -1.8  x  10  per  datum  or  roughly  -10  per  uncprrelated 

noise  value.  Owing  to  the  clutter-edge  effect,  it  is  difficult  to 

estimate  the  critical  sample  size  in  clutter  from  Fig.  4.  However, 

if  £  has  the  same  value  in  clutter  as  in  uniform  sky,  then  N  ~  150 

c 

data  (6  mR)  in  clutter,  which  agrees  well  enough  with  the  appearance 
of  the  curves  in  Fig.  4.  The  critical  sample  size  is  always  greater 
in  uniform  sky  than  in  clutt  r,  provided  £  is  the  same  in  the  two 
regions,  because  |a|  is  always  greater  in  uniform  sky  than  in 
clutter. 
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B.  Threshold  Parameter,  Critical  Sample  Size,  and 
Clutter-edge  Effect 

The  sampling  s.d.s  of  the  uniform  sky  and  clutter  populations  of 
Figs.  3E  and  3F  are  replotted  against  frame  size  in  Fig.  4A. 

Obviously  in  both  cases  the  sampling  s;d.  decreases  at  a  decreasing 
rate  with  increasing  frame  size.  This  is  not  surprising;  the  law  of 
diminishing  returns  with  increasing  sample  size  is  well  known.  More 
important  is  the  relation  between  the  sampling  s.d.  and  the  threshold 
parameter. 

To  study  this  relation  we  first  calculated  threshold  parameter 

ratios  x  /<|x  |>  from  the  two  data  sets  (each  containing  114  segments) 
w  r 

used  in  the  previous  article,^  not  from  the  restricted  sets  of 
Figs.  3E  and  3F.  (The  frame  sizes  employed  were  2n  ■  36,  80,  144, 

240,  and  340  data,  the  limit  of  340  being  necessary  to  prevent 
excessive  loss  of  ratios.  A  pair  of  data  added  to  the  frame  reduces 
the  number  of  ratios  calculable  from  a  segment  by  two.  The  smallest 
of  the  228  segments  each  contained  I486'  data.  A  frame  size  of  340 
data  is  23  percent  of  a  smallest  segment,  which  means  that  23  percent 
.of  the  data  in  those  segments  could  not  be  used  as  window  values  x^. 

If  larger  frame  sizes  had  been  used,  the  results  would  be  unreliable 
because  of  the  additional  ratios  lost.)  Next  we  arranged  the  ratios 
in  statistical  frequency  tables  from  which  we  determined  threshold 
parameters  that  allow  one  false  alarm  (f.a.)  in  all  the  data  or  one 
in  the  uniform  sky;  these  parameters  are  plotted  in  Fig.  4B.  The  one 


stationary  noise.  If  the  noise  is  correlated  and  non-stationary ,  B 
is  determined  by  these  two  properties  jointly,  and  a  numerical 
analysis  would  require  calculations  with  synthetic  noise  having  a 
known  type  and  amount  of  non-stationarity .  The  effort  did  not  seem 
worthwhile.  In  this  case  we  simply  regard  equation  (2)  as  an 
empirical  relation  useful  for  interpreting  and  predicting  the  effects 
of  changing  threshold-sensor  size. 

This  treatment  of  non-stationarity  needs  clarification,  lest  it 
seem  flawed.  Non-stationary  noise  has  no  definite  mean  absolute 
deviation,  so  the  sampling  variance  of  that  statistic  cannot  approach 
zero  and  may  not  decrease  raonotonically  with  increasing  sample  size. 
True,  for  an  infinite  data  set,  but  real  data  sets  are  finite.  A 
finite  data  set  does  have  a  definite  mean  absolute  deviation  whose 
value  is  obtained  without  error  when  the  sample  includes  all  data  of 
the  set.  An  infinite,  non-stationary  data  set  is  equivalent  to  an 
infinite  ensemble  of  finite  data  sets,  each  having  different 
statistics.  One  always  works  with  a  few  members  of  the  ensemble  at 
most.  Conclusions  thus  reached  will  be  useful  if  the  data  fairly 
represent  the  conditions . the  surveillance  set  encounters. 


log( o / o )  -  B  +  a  log(/N).  <2) 

The  line  through  the  uniform-sky  results  is  adequately  fit  by  a  « 
-0.909  and  8  *  0.0819,  that  through  the  clutter  results  by  a  ■  -0.530 
and  6  *  -0.0433. 

The  fact  that  a  >  -1  for  the  uniform-sky  results  requires 
comment.  The  uniform-sky  noise  may  not  be  strictly  stationary 
because  small  gain  differences  (•<  8%)  between  segments  were  not 
corrected.  (Neither  were  they  corrected  for  the  simulations  here  or 
in  the  previous  article.^)  On  the  other  hand,  non-stationarity 
between  segments  should  increase  the  ranges  of  the  sampling  a.d.s 
primarily,  whereas  non-stationarity  within  the  segments  should  be  the 
chief  cause  of  increases  in  the  slopes.  The  small  differences 
between  -1  and  a  for  the  uniform-sky  results  may,  therefore,  be  due 
to  ordinary  experimental  uncertainty  like  that  in  Figs.  3A  and  3B, 
not  to  gain  differences. 

If  the  noise  is  stationary,  B  can  be  related  to  the  fraction  f 

of  the  data  which  are  effectively  independent.  The  lines  through  the 

sampling  s.d.s  of  the  correlated  Gaussian  noise  have  a  *■  -1  and  B  “ 

0.188.  From  the  latter  we  .find  f  ■  0.42.  Calculations  of  product- 

moment  correlation  coefficients  show  that  the  first  and  last  of  five 

sequential  magnitudes  in  the  LMS-filtered  synthetic  noise  are 

uncorrelated,  giving  f  m  0.40.  The  theoretical  correlation  length 

1 3 

of  LMS-filtered  white  noise  gives  -  0.5  approximately. 

Clearly  the  suggested  interpretation  of  B  is  sound  for  correlated 


to  obtain  Gaussian  noise  with  the  same  correlation  as  the  LMS- 
filtered  uniform-sky  noise.  The  correlated  synthetic  noise  was  then 
processed  in  the  same  way  as  the  uncorrelated  noise;  results  are 
plotted  in  Figs.  3C  and  3D.  Correlation  and  non-random  sampling 
together  reduce  the  number  of  independent  values  in  the  frame,  so 
that  the  values  and,  usually,  the  ranges  of  the  sampling  s.d.s  are 
increased.  For  reasons  already  given  the  sampling  s.d.s  at  a  frame 
sire  of  eight  data  are  closer  to*  the  diagonals  than  those  at  larger 
frame  sizes,  which  alone  are  of  practical  importance.  At  the  larger 
frame  sizes  correlation  and  non-random  sampling  simply  shift  the 
reference  line  rightward  parallel  to  itself,  as  illustrated  by  the 
lines  with  slope  -1  drawn  through  the  points  at  frame  sizes  above 
eight  data,. 

Finally,  sampling  variances  were  calculated  for  ten  segments  of 
uniform-sky  noise  (each  1900  values  long)  and  ten  of  clutter  (each 
3200- values  long).  Comparing  these  results  -  Figs.  3E  and  3F  -  with 
those  for  correlated  synthetic  noise,  we  see  that  non-random  sampling 
and  non-stationarity  increase  the  ranges  of  the  sampling  s.d.s  and 
the  slopes  of  the  lines  through  them,  most  obviously  in  the  case  of 
the  clutter.  The  sampling  s.d.s  of  both  the  uniform-sky  noise  and 
the  LMS-filtered  clutter  can  be  described  empirically  at  the  frame 
sizes  of  Interest  simply  by  extending  equation  (1)  to  include  the 
effects  of  non-random  sampling,  correlation,  and  non-stationarity. 

The  general  equation  for  any  line  in  Fig.  3  is 


correlated  and  non-statlonary .  (Correlation  i6  measured  in  this  work 
by  the  product-moment  correlation  coefficient.**)  To  relate  the 
theorem  to  the  problem  in  hand,  we  made  computations  in  which  we 
first  approximated  the  conditions  of  the  theorem  and  then  relaxed 
them  by  steps  to  arrive  at  the  actual  conditions. 

Nine  segments  of  uncorrelated,  Gaussian,  pseudo  random  noise 
were  generated  by  the  method  described  in  reference  12.  This  was 
done  for  two  segment  lengths,  1900  and  3200  values,  which  are  the 
average  lengths  of  the  segments  of  recorded  uniform-sky  noise  and 
clutter  used  in  the  previous  article.*  Sampling  variances  were 
calculated  by  passing  the  segments  through  the  adaptive  threshold 
algorithm  using  ten  frame  sizes  from  8  to  640  values.  The  results 
appear  in  Figs.  3A  and  3B,  where  the  dots  indicate  the  sampling 
s.d.s  of  the  whole  data  set,  and  the  bars  span  the  ranges  of  the 
sampling  s.d.s  of  the  segments.  (The  variances  are  calculated  about 
H  ■  <|x|>,  the  average  absolute  deviation  of  all  data  in  a  set. 
Consequently,  a  population  sampling  variance  is  the  average  of  the 
associated  segment  sampling  variances.)  The  sampling  s.d.s  for  the 
population  are  on  or  very  near  the  diagonals  showing  that  sampling  by 
an  adaptive  threshold  sensor  is  equivalent  to  random  sampling  if  the 
noise  is  stationary  and  uncorrelated.  The  ranges  of  the  sampling 
s.d.s  for  the  segments  result  from  chance  differences  between  the 
segments.  They  are  larger  for  the  shorter  segments  as  expected. 

Next  the  synthetic  noise  was  passed  through  the  1-D  LMS  filter 


A.  Numerical  Analysis  of  Sampling 


Since  dispersion  is  measured  by  the  mean  absolute  deviation,  the 
sensor's  sampling  performance  can  be  related  to  the  theorem  on  the 
sampling  distribution  of  the  meanj  which  states  in  part:^  If  random 
samples  of  2n  ■  N  values  each  are  taken  from  a  population  (data  set) 


of  more  than  20N  values,  the  relation  between  the  population  variance 


o  and  the  sampling  s.d.  of  the  mean  is  closely  approximated  by 


o *  o^/N  or  log (a^/a)  -  -  logCVN).  (1) 

n  h 


The  sampled  population  consists  of  the  absolute  deviations  |x  -  p|, 
which  are  the  noise  magnitudes  |xj  since  the  LMS-filtered  data  have 
zero  means  locally  and  globally. 

The  logarithmic  form  of  equation  (1)  is  represented  by  the 
diagonals  in  Fig.  2.  The  sampling  s.d.  is  divided  by  the  population 
s.d.  to  facilitate  comparison  of  results  from  data  sets  with 
different  variances.  It  is  useful  to  think  of  a  one-sided  threshold 
sensor  and  consider  the  limit  of  equation  (1)  when  N  =  1.  In  that 
.case  h  -  |x,.|  so  that  <J„  -  a  by  definition,  and  the  equation  is 

f  n 

an  identity.  Thus,  we  expect  equation  (1)  to  be  nearly  satisfied 
when  the  sample  size  is  small,  regardless  of  other  conditions. 

Except  for  N  ■  1 ,  however,  the  theorem  is  not  applicable  to  the 
adaptive  threshold  sensor  and  the  recorded  background  noise  because 
the  sensor  does  not  take  random  samples,  and  the  recorded  noise  is 


increasing  frame  size,  but  the  threshold  parameter  may  decrease  or 
increase,  and  signal  detection  may  improve  or  deteriorate  depending 
on  other  factors.  To  transform  this  qualitative  description  into  a 
semi-quantitative  one,  we  analyzed  the  threshold  sensor's  sampling 
performance,  related  that  to  the  threshold  parameter  required  for  a 
given  false  alarm  rate,  and  finally  linked  the  threshold  parameter 
signal  detection. 


the  sensor  of  Fig-  1  cannot  operate  statistically;  its  behavior  in 
that  case  is  described  in  the  appendix. 


(Here  a  digression  is  necessary  to  explain  an  important  detail. 

With  uniform-sky  noise,  which  is  stationary,  the  local  values  T)^  of 

the  mean  absolute  deviation  are  the  same  as  its  global  value  n  .  Th 

£ 

frame  averages  n  estimate  both  tj.  and  f|  .  The  variance  of  fj  about 

1  S 

is  due  only  to  chance  differences  between-  samples.  With  a  finite 

set  of  non-stationary  clutter  data,  the  local  values  differ  from 

the  global  value  n  *  <|x|>  for  the  whole  data  set.  The  frame 
g 

averages  still  estimate  both  n,  and  n  ,  and  chance  differences 

1  g 

between  samples  still  contribute  to  the  variance  of  n  about  n  *  but 

g 

the  variance  is  Increased  by  the  non-stationarity . ) 

This  summary  of  the  adaptive  threshold  sensor's  operation 
implies  certain  effects  which  occur  as  the  frame  size  is  increased. 
When  the  frame  is  small  some  estimates  of  r|  ■  <|x|>  are  much  too 
small,  others  much  too  great.  Large  fj  degrade  signal  detection 
locally  by  causing  the  local  threshold  Kn  to  be  larger  than 
necessary  to  prevent  a  false  alarm.  Small  n  degrade  signal 
detection  globally  since  the  threshold  parameter,  K  *  limit  of 
(x^/n) ,  is  the  same  throughout  the  scene  and  must  be  large  to 
prevent  a  false  alarm  where  rj  is  small.  As  the  frame  size  is 
increased  the  estimates  of  n  become  more  accurate,  the  sampling 
variance  and  the  threshold  parameter  decrease,  and  signal  detection 
improves,  provided  the  noise  is  stationary.  If  the  noise  is 
non-stationary  the  sampling  variance  may  still  decrease  with 
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shows  cloud-edge  clutter  at  the  outputs  of  the  LMS  filter  and  the 
wideband  filter  used  for  data  collection.  In  this  work  the  optical 
sensor  is  the  one  used  for  data  collection,  but  the  wideband  filter 
is  virtually  eliminated  by  passing  the  data  through  the  LMS  filter  to 
obtain  a  low  value  of  <|x|>  throughout  the  scene.  (This  is 
permissible  because  the  LHS  passband  lies  well  within  that  of  the 
wideband  filter.)  With  the  optical  sensor  and  electronic  filter 
specified,  the  quality  of  signal-  detection  depends  mainly  on  the 
value  of  K  required  for  the  designated  false  alarm  rate. 

The  mean  p  of  the  filtered  noise  determines  the  mode  of  the 

sensor's  operation  as  follows:  If  the  mean  is  always  zero  over  the 

sensor's  reach  ("locally"),  it  is  also  zero  over  all  greater  extents 

("globally").  The  mean  of  uniform-sky  noise  is  zero  locally  and 

globally,  as  is  the  mean  of  LMS-filtered  clutter  (Fig.  2B) .  In  such 

cases  the  mean  magnitude  <|x|>  and  mean  absolute  deviation  <[x  -  p|> 

are  identical  and  are  estimated  locally  and  globally  by  the  frame 

averages  fi  *•  <  |  x^.  |  >  -  The  estimate  q  and  the  local  threshold 

(p  +  Kq)  -  ICq  vary  as  the  filter  output  passes  through  the  sensor. 

The  distribution  and  variance  of  q  about  q  =  <|x|>  are  known  as  the 

sampling  distribution  and  sampling  variance;  the  standard  deviation 

9 

(s.d.)  is  known  as  the  standard  error  of  the  mean,  but  we  shall  call 
it  the  sampling  s.d..  With  noise  whose  mean  is  zero  locally  and 
globally,  the  sensor  of  Fig.  I  operates  "statistically",  i.e.,  sets 
the  threshold  at  T  *■  (p  +  Kfi)  and  makes  a  statistical  test  of  the 
window  value's  signal  content.  If  the  noise  mean  is  not  zero  locally. 
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Because  of  difficulties  in  estimating  the  local  probability 
distribution  and  its  parameters,  no  adaptive  threshold  sensor 
operates  ideally.  The  simpler,  practical,  1-D  sensor  studied  here  is 
represented  in  Fig.  1  with  discrete  values  of  the  filter  output 
passing  stepwise  through  it.  At  feach  step  this  sensor  averages  the 
magnitudes  of  2n  values  which  lie  equally  about  the  single  value 
whose  signal  content  is  to  be  tested.  The  local  threshold  is  a 
constant  multiple  K  of  the  average  <|x^|>.  The  part  of  the  sensor 
containing  the  tested  value  will  be  called  the  signal  window,  and  the 
part  containing  the  2n  averaged  magnitudes  will  be  called  the  frame. 
The  part  called  the  shutter  blocks  values  from  the  wings  of  a 
point-source  signal  centered  on  the  window  and  prevents  them  from 
affecting  the  threshold.  Since  the  threshold  parameter  K  is 
constant,  the  false  alarm  probability  varies  with  the  local 
probability  distribution  as  well  as  at  clutter  edges.  Nonetheless,  K 
can  be  made  large  enough  to  provide  any  desired  false  alarm  rate  in 
the  whole  scene. 

The  mean  magnitude  <|x|>  of  the  filtered  noise  determines  the 
quality  of  this  sensor's  operation.  The  lower  the  average  threshold 
K<|x|.>,  the  better  the  signal  detection  is  on  average.  Ordinarily 
< | x | >  affects  signal  detection  more  than  K  because  noise  magnitude 
can  vary  as  much  as  20  dB  with  different  optical  sensors  and 
electronic  filters;  the  range  of  the  threshold  parameter  is  smaller 
unless  there  is  an  extreme  clutter-edge  effect.  Differences  in  noise 
magnitude  with  different  filters  are  illustrated  in  Fig.  2  which 


II.  Results  and  Discussion 


An  adaptive  threshold  sensor  rejects  noise  and  detects  signals 
through  a  statistical  test  of  the  null  hypothesis,  filter  output  * 
noise,  and  the  alternative,  filter  output  -  (signal  +  noise). 

Ideally  the  sensor  does  this  by:  taking  a  sample  of  N  values  from  the 
local  noise;  estimating  the  local  mean,  dispersion,  and  probability 
distribution  from  the  sample;  setting  a  threshold  t  from  the 
estimates;  taking  a  sample  of  one  value  x^  from  the  local  noise  not 
previously  sampled;  rejecting  the  null  hypothesis  if  x  >  £.  A 
threshold  crossing  is  either  a  signal  or  a  false  rejection  of  the 

g 

null  hypothesis  (false  alarm).  Further  discriminants  are  needed  to 
distinguish  the  two  under  operating  conditions,  but  not  in  a 
simulation  because  it  is  knovm  whether  a  simulated  signal  was  added 
to  the  background  data.  The  dispersion  is  usually  measured  by  the 
standard  deviation  a  or  the  mean  absolute  deviation  n.  The  local 
threshold  is  the  sum  of  the  local  mean  p  and  a  multiple  of  the  local 
dispersion,  e.g.,  t  -  (p  +  kn) .  The  multiple  is  chosen  to  provide  a 
desired  probability  of  false  alarm  and  so  depends  on  the  local 
probability  distribution.  Nevertheless,  the  false  alarm  probability 
will  vary  at  the  "clutter  edges"  described  in  section  11(B), ^  and 
these  can  occur  if  either  the  distribution  or  its  parameters  alone 
vary.  Thus,  constant  false  alarm  probability  is  unattainable  in  a 
cluttered  scene  even  with  a  statistically  ideal  adaptive  threshold 


sensor . 


LMS-filtered  noise.  This  leads  to  some  useful  experimental 
generalizations  and  to  a  principle  for  determining  by  data  analysis 
whether  a  2-D  threshold  sensor  may  be  advantageous.  The  methods 
employed  are  applicable  wherever  recorded  background  data  are 
available  and  should  be  useful  for  developing  threshold  sensors  to 
operate  under  conditions  different  from  the  ones  in  this  simulation. 
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mathematically  tractable.  Our  simulation  incorporates  the 
theoretical  ideas  but  is  applicable  to  real  cases  because  recorded 
natural  background  noise  is  used  instead  of  clutter  models  and  no 
simplifying  assumptions  need  be  made  -  the  algorithms  can  be 
programmed  to  operate  on  the  noise  data  as  the  filter  and  threshold 
sensor  would  operate  on  the  sampled  output  of  the  optical  sensor. 

(For  convenience  the  sensor  algorithm  actually  was  simplified  by 
ignoring  off-peak  detection  of  signals.  The  resulting  error  is 
negligible,  as  determined  by  computation.  See  the  last  paragraph  of 
section  11(C)  for  details.) 

Because  recorded  background  noise  and  exact  algorithms  are  used, 
the  simulation  provides  inputs,  intermediate  quantities,  and  outputs 
of  the  simulated  devices  under  actual  operating  conditions.  By 
studying  the  properties  and  interrelations  of  the  inputs, 
intermediates,  and  outputs  while  varying  the  device  parameters,  one 
gains  both  insight  into  the  operation  of  the  devices  and  quantitative 
design  information.  In  this  connection,  it  should  be  noted  that  the 
simulations  of  this  article  and  the  previous  one^  emphasize 
statistical  analysis  of  the  data  at  all  stages;  they  do  not  simply 
mimic  device  behavior.  In  .fact,  a  major  purpose  of  this  work  is  to 
develop  techniques  of  simulation  which  incorporate  statistical 
analysis  for  depign  purposes,  and  to  demonstrate  the  value  of  these 
techniques  by  example.  With  that  in  mind,  we  present  a  numerical 
study  of  the  threshold  sensor's  sampling  performance,  taking  account 
of  non-random  sampling,  correlation,  and  non-stationarity  of  the 


The  plot  of  sampling  s.d.  in  clutter  suggests  that  the  threshold 
parameter  for  one  false  alarm  in  clutter  should  decrease  slowly 
beyond  *  150  data.  In  fact,  it  passes  through  a  broad,  shallow 
minimum  because  of  the  following  clutter-edge  effect,  which  appears 
to  begin  at  a  frame  size  between  8  and  10  mR  (200  and  250  data).  As 
the  size  of  the  threshold  sensor  is  increased,  the  frame  includes 
larger  and  larger  proportions  of  the  benign  clutter  which  flanks 
severe  clutter  like  that  from  1 58  to  165  mR  in  Fig.  2B.  This  lowers 
the  estimates  fj  *  <|x^|>  and  raises  the  ratios  x  /rj  in  the 
vicinity  of  the  severe  clutter,  thus  increasing  the  threshold 
parameter  needed  to  prevent  a  false  alarm.  Ordinarily  the  increase 
will  be  no  more  than  a  few  decibels;  only  an  extreme  clutter-edge 
effect  like  the  hypothetical  one  described  in  section  V(B)  of  the 
previous  article^  could  change  the  threshold  parameter  by  10-20  dB. 
Partly  for  this  reason,  the  threshold  sensor  usually  affects  signal 
detection  less  than  the  optical  sensor  and  electronic  filter,  as 
noted  earlier.  Here  the  increase  of  the  threshold  parameter  -  2.3% 
between  frame  sizes  of  240  and  340  data  -  is  within  the  experimental 
uncertainty.  However,  the  occurrence  of  a  clutter-edge  effect  fits 
ao  well  into  the  pattern  of  the  results  that  we  do  not  doubt  its 
actuality.  With  that  granted  the  dashed  part  of  the  top  curve  in 
Fig.  4B  indicates  how  the  threshold  parameter  might  behave  at  larger 
frame  sizet.. 

To  verify  or  refute  this  behavior  would  require  computations 
with  unsegraented  data.  Many  frame  size  effects  can  be  studied 
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without  separating  uniform-sky  noise  from  clutter,  and  if  the  data 
were  not  segmented,  much  larger  frame  sizes  could  be  used  before  the 
loss  of  window  values  mentioned  above  became  a  problem.  This  would 
be  advantageous  for  frame  size  studies,  but  would  hide  differences 
between  performance  in  uniform  sky  and  clutter.  We  chose  to 
emphasize  the  latter,  realizing  that  frame  size  studies  with 
separated  data  are  valid  provided  the  clutter  segments  do  not  begin 
or  end  at  clutter  edges,  a  condition  satisfied  here. 

(The  clutter-edge  effect  has  also  been  called  the  cloud-edge 
effect,  but  bright  cloud  edges  may  not  be  its  only  cause.  The  term, 
clutter-edge,  is  preferable  since  it  accurately  describes  the  origin 
of  the  effect  in  the  waveform  but  is  noncommittal  about  its  origin  in 


C.  Signal  Detection,  Optimum  Sample  Size,  and  Beat  Sample  Size 

Frequencies  of  detecting  simulated  signals  were  calculated  using 
the  threshold  parameters  in  Fig.  4B  and  the  two  data  sets  of  114 
segments  each  from  uniform  sky  and  from  clutter.  The  computed  values 
are  frequencies  of  detecting  simulated  signals  at  their  peaks.  There 
is  an  additional  small  chance  for  signals  missed  at  their  peaks  to  be 
detected  elsewhere  on  the  waveform.  Ignoring  off-peak  detection 
causes  <  2.5%  error  in  detection  frequencies  >  0.5,  as  determined  by 
calculation.  Specifically,  if  the  frequency  of  detection  at  the  peak 
is  0.5000,  the  actual  detection  frequency  including  off-peak 
detection  is  <  0.5125.  Frequency  of  off-peak  detection  decreases 
with  increasing  frequency  of  detection  at  the  peak,  since  only  one 
detection  of  a  signal  is  independent  and,  therefore,  off-peak 
threshold  crossings  are  not  counted  when  the  signal  is  detected  at 
its  peak.  If  off-peak  detection  is  ignored,  the  adaptive  threshold 
algorithm  is  much  simpler  and  faster,  yet  the  resulting  error  is 
negligible.  Accordingly,  this  approximation  was  used  here  and  in  the 
previous  work.^ 

•Frequencies  of  detection  at  the  peak  are  plotted  in  Fig.  5 
except  those  for  the  frame  size  of  1^4  data,  which  are  in  the 
previous  article^  and  are  hardly  distinguishable  from  the  results  for 
the  frame  size  of  340  data.  From  these  curves  we  obtained  data  for 
Fig.  6,  where  simulated  signal  strengths  required  for  detection 
frequencies  of  0.50,  0.70,  0.90,  and  0.95  are  plotted  against  frame 


size.  The  ordinates  of  curves  in  Fig.  6  will  be  called  equally 
detected  signal  strengths. 

The  curves  of  equally  detected  signal  strengths  are  much  like 
the  top  threshold  parameter  curve  in  Fig.  4B,  indicating  that 
threshold-sensor  size  influences  signal  detection  chiefly  through  its 
effect  on  the  threshold  parameter.  The  minima  of  the  curves  occur, 
regardless  of  detection  frequency,  at  a  frame  size  of  about  10  mR 
(250  data),  the  same  as  the  minimum  of  the  threshold  parameter  curve. 
This  is  the  optimum  sample  size  for  signal  detection.  Increasing  the 
frame  size  to  14  or  15  mR  increases  equally  detected  signal  strengths 
about  0.2  dB,  the  same  as  the  causative  increase  of  the  threshold 
parameter.  Because  of  the  scales  employed  in  Fig.  6,  signal 
detection  appears  to  deteriorate  quickly  below  the  optimum  size,  but 
actually  it  does  not.  Using  frames  of  144  or  36  data  rather  than 
the  optimum  size  reduces  signal  detection  by  only  0.2  dB  and  2.4  dB 
respectively.  Thus,  the  best  sample  size  here  is  the  critical  sample 
size  in  clutter  since  use  of  a  larger  frame  risks  a  clutter-edge 
effect  but  improves  signal  detection  very  little  in  either  clear  sky 
or  clutter.  The  frame  must  not  be  too  much  less  than  the  optimum 
size,-  however,  because  the  curves  in  Figs.  4B  and  6  are  rising 
rapidly  at  36  data,  one-seventh  the  optimum  size. 


D.  Reality  of  the  Simulation  and  Limitations  of  the  Results 


These  subjects  are  discussed  at  length  in  the  previous  article 
(sections  I,  V(B),  VI  and  appendix  A),^  but  it  may  be  worthwhile  to 
re-emphasize  two  points  here. 

The  first  concerns  the  "reality"  of  the  simulation.  There  are 
two  realities,  that  of  the  signal  processing  and  that  of  the  scene 
and  signal  sources.  As  regards  -signal  processing,  the  simulation  and 
reality  are  alike  since  the  algorithms  effectively  operate  on  the 
noise  data  as  the  simulated  devices  would  operate  on  the  output  of 
the  optical  sensor.  However,  the  relation  of  simulated  signals  to 
sources  in  the  scene  departs  from  reality.  The  simulated  signals 
should  correspond  to  sources  of  constant  radiance  placed  at  each 
sampled  position  of  the  scene,  with  account  taken  of  the  background 
covered  by  the  sources  (reference  1,  appendix  A).  That  was  impossible 
in  this  work  because  the  data  could  not  be  radiometrically  calibrated. 
Consequently,  the  simulated  signals  used  here  correspond  to  real 
sources  of  constant  "excess  radiance"  placed  at  each  sampled  position 
of  the  scene.  ("Excess  radiance"  is  radiance  above  that  of  the 
background  covered  by  the  source,  all  radiances  being  measured  at  the 
source  plane.)  In  clear  sky  where  the  background  radiance  is 
uniform,  signals  from  the  two  types  of  sources  are  the  same,  but  they 
are  different  in  clutter  where  the  background  radiance  varies.  To 
recapitulate,  the  simulated  signals  added  to  the  uniform-sky  noise 
and  to  the  clutter  correspond,  respectively,  to  real  sources  of 


constant  radiance  and  to  real  sources  of  constant  excess  radiance.' 

For  that  reason,  the  two  types  of  simulated  signals  are  called, 
respectively,  signals  and  pseudo  signals.  Use  of  pseudo  signals 
makes  the  computed  frequencies  of  detection  too  large  in  clutter, 
i.e.,  the  detection  frequencies  in  Fig.  5B  actually  require  larger 
signal  strengths  than  those  shown.  Thus  the  simulated  performance  is 
better  in  clutter  than  the  actual  performance  by  an  unknown  and,  with 
the  given  data,  an  unknowable  amount.  This  should  not  significantly 
•affect  comparison  of  signal  processing  options,  however,  because  all 
options  are  tested  against  simulated  signals  with  the  same  connection 
to  real  sources. 

The  second  point  for  re-emphasis  concerns  limitations  of  the 
results.  Since  this  simulation  is  a  form  of  data  analysis,  the 
results  are  restricted  by  the  data  employed,  specifically,  by  all  the 
natural  and  man-made  conditions  existing  when  the  data  were 
collected.  Relative  to  the  present  simulation  the  most  important 
experimental  conditions  are:  (Natural)  time  and  place  of 
measurements,  meteorologic  conditions,  'type  and  range  of  clouds; 
(Man-made)  line  of  sight,  optical  resolution,  spectral  passband, 
dwell  time  of  the  point-source  image  on  the  detector,  point-source 
output  of  the  detector,  electronic  passband  and  dynamic  range,  NEI, 
sampling  rate,  and  sampling  time.  Section  II  of  the  previous 
article^  and  the  appendix  to  this  article  give  all  available 
information  about  these  conditions;  any  such  information  not  found  in 
either  place  is  not  available. 


Despite  restrictions  imposed  by  experimental  conditions,  data 
analysis  often  leads  to  inferences  of  more  general  validity,  and  that 
is  true  here.  In  the  next  section,  the  more  general  and  the  more 
restricted  results  are  clearly  disjoined,  and  conditions  for  the 
validity  of  each  are  plainly  stated.  The  conditions  are  inseparable 
from  the  results  and  conclusions;  the  latter  must  not  be  read, 
recalled,  or  used  without  the  former. 


III.  Summary  and  Conclusions 


We  first  summarize  results  and  state  conclusions  which  are  valid 
in  general  for  the  type  of  adaptive  threshold  sensor  in  Fig.  1,  when 
it  follows  a  spatial  filter  whose' mean  output  is  zero  over  extents  2 
the  sensor's  reach.  Barring  an  extreme  clutter-edge  effect,  the 
quality  of  signal  detection  depends  primarily  on  the  average 
magnitude  of  the  filtered  noise  -(set  by  the  optical  sensor  and 
electronic  filter)  and  secondarily  on  the  threshold  parameter.  The 
threshold  parameter  required  for  a  given. false  alarm  rate  is 
determined  by  the  threshold  sensor's  sampling  variance  and  by  a 
clutter-edge  effect.  Ordinarily  these  factors  change  the  threshold 
parameter  and  signal  detection  no  more  than  a  few  decibels  within  the 
range  of  practical  threshold  sensor  sizes.  Only  an  extreme 
clutter-edge  effect  could  result  in  changes  greater  than  10  dB. 
Sampling  variance  decreases  at  a  decreasing  rate  with  increasing 
threshold-sensor  (sample)  size,.  The  decrease  is  more  rapid  when  the 
noise  is  stationary  than  when  it  is  non-stationary .  There  are 
loosely  delimited  critical  sample  sizes  above  which  sampling  variance 
decreases  very  slowly.  The  critical  sample  size  for  stationary  noise 
is  about  100  independent  data  at  both  radar  and  IR  wavelengths.  The 
critical  sample  size  is  larger  for  clear-sky  (stationary)  noise  than 
for  clutter  (non-stationary)  noise.  If  the  scene  is  uncluttered,  the 
threshold  parameter  and  equally  detected  signal  strength  behave  as 
the  sampling  variance  does.  If  the  scene  is  cluttered,  they  behave 


likewise,  provided  there  is  no  clutter-edge  effect.  If  there  is  one, 
the  threshold  parameter  for  a  given  false  alarm  rate  in  the  whole 
scene  passes  through  a  minimum  with  increasing  sensor  size,  causing 
equally  detected  signal  strength  to  do  the  same  in  both  the  clear  and 
cluttered  parts  of  the  scene,  as  in  Figs.  4B  and  6. 

The  critical  sample  size  is  as  large  as  the  adaptive  threshold 
sensor  of  Fig.  1  need  be,  since  making  it  larger  hardly  improves 
performance.  On  the  other  hand,* the  sensor  cannot  be  larger  than  the 
onset  of  a  clutter-edge  effect  without  loss  of  performance.  Thus  if 
a  clutter-edge  effect  occurs,  it  may  be  impossible  to  expand  the 
sensor  enough  in  the  scan  direction  to  reach  the  critical  sample  size 
in  clear  sky  or  even  in  clutter.  With  a  column  of  detectors  normal 
to  the  scan,  however,  sampling  in  two  directions  is  possible.  This 
allows  the  sample  size  to  be  increased  without  expanding  the  sensor's 
length,  albeit  at  the  risk  of  clutter-edge  effects  in  all  directions 
within  the  field  of  view.  The  question  then  is  whether  1-D  and  2-D 
sensors  perform  equally  or  whether  one  type  is  better. 

The  following  principle  will  be  useful  for  investigating  that 
question:  If  the  onset  of  a  clutter-edge  effect  with  a  1-D  sensor  is 

less  than  the  critical  sample  size  in  either  clutter  or  clear  sky,  a 
2-D  sensor  may  be  (but  is  not  guaranteed  to  be)  advantageous, 
otherwise  it  will  not  be.  Note  that  if  a  1-D  sensor  fails  to  achieve 
the  critical  sample  size  in  clutter,  a  2-D  sensor  may  do  so,  and  thu9 
be  advantageous,  even  if  it  too  fails  to  reach  the  critical  sample 
size  in  clear  sky.  The  intended  use  of  this  principle  is  evaluation 


of  threshold  sensor  performance  by  statistical  analysis  of  background 
data.  To  apply  the  principle,  one  must  know  the  critical  sample 
sizes  and  the  onset  of  the  clutter-edge  effect.  The  critical  sample 
sizes  can  be  estimated  from  the  sampling  variance,  the  threshold 
parameter,  or  signal  detection.  The  sampling  variance  gives  the 
least  accurate  estimates  but  is  by  far  the  easiest  method  to  use. 

The  clutter-edge  effect  and  its  onset  can  be  found  only  from  the 
threshold  parameter  or  signal  detection.  The  principle  i  limited  to 
cases  where  the  sensor  makes  a  statistical  test  for  signals,  i.e., 
cases  where  the  noise  mean  is  zero  over  extents  2  the  sensor's  reach. 
This,  requires  the  filter  preceding  the  sensor  to  have  a  high  cut-on 
frequency.  The  principle  is  not  applicable  with  a  low-pass  filter  or 
a  bandpass  filter  having  a  low  cut-on  frequency.  The  clutter  outputs 
of  these  filters  simply  cause  the  sensor  to  turn  off  the  surveillance 
set,  as  explained  in  the  appendix. 

(Some  hypothetical  cases  can  be  described  to  exemplify  the 
principle.  If  a  1-D  and  a  2-D  sensor  both  achieve  the  critical 
sample  size  in  clear  sky  without  a  clutter-edge  effect,  they  perform 
equally.  That  is  true  because  they  estimate  the  local  noise 
statistics  in  the  same  way  ‘from  samples  of  the  same  size  and  test  the 
same  statistical  hypotheses.  A  2-D  sensor  may  perform  better  than  a 
1-D  sensor  if  the  clutter  is  spatially  homogeneous  and  the  1-D  sensor 
incurs  a  clutter-edge  effect  before  reaching  both  critical  sample 
sizes.  In  this  case  a  2-D  sensor  may  achieve  one  or  both  critical 
sample  sizes  without  a  clutter-edge  effect.  It  is  necessary  to  be 


more  precise  about  the  meaning  of  "homogeneous",  since  clutter  is 
intrinsically  heterogeneous.  in  the  present  context,  clutter  is 
homogeneous  if  its  statistical  properties  change  slowly  in  all 
directions  over  extents  £  the  sensor's  reach.  If  the  clutter  is 
spatially  inhomogeneous,  a  2-D  sensor  might  be  inferior  to  a  1-D 
sensor.  This  could  occur  if  the  noise  statistics  changed  slowly  in 
the  scan  direction  but  rapidly  in  another  direction,  resulting  in  a 
clutter-edge  effect  with  the  2-D.  but  not  the  1-D  sensor.  A  bright 
cloud  edge  parallel  to  the  _^can  might  lead  to  that  outcome.  The 
principle  applies  in  all  these  cases  because  it  involves  only 
invariable  relations  between  statistical  properties  of  the  noise  and 
sensors.  The  principle  is  inapplicable  (e.g.,  with  low-pass  clutter) 
only  if  the  noise  and  sensors  do  not  have  these  properties.  However, 
the  intended  use  of  the  principle  is  analysis  of  data,  not 
hypothetical  cases.  The  principle  tells  one  how  to  determine  whether 
1-D  and  2-D  threshold  sensors  perform  differently,  but  the  actual 
determination  awaits  analysis  of  data.) 

We  next  state  less  general  results'  and  conclusions.  These  are 

valid  only  for  a  1-D  adaptive  threshold  sensor  following  a  1-D 

LMS  filter  when  the  filter  .input  is  from  an  optical  sensor  having  an 
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NEI  2:  1.5  x  10  W/ cm  and  a  spatial  resolution  better  than 
0,15  x  0.36  mR.  Further  restrictions  are  imposed  by  all  the  natural 
and  man-made  experimental  conditions  existing  when  the  data  were 
collected,  as  discussed  above  in  section  11(D)  and  summa.ized  in  the 


appendix . 


Under  the  stipulated  conditions,  the  standard  critical  sample 
size  in  clear  sky  (100  independent  values)  corresponds  to  250  data 
extending  over  10  mR;  the  critical  sample  size  in  clutter  is  6  mR  or 
150  data.  A  clutter-edge  effect  with  an  onset  between  8  and  10  mR 
results  in  an  optimum  threshold-sfensor  size  of  10  mR,  but  increasing 
the  size  to  14-15  mR  increases  the  threshold  parameter  and  equally 
detected  signal  strengths  only  0.2  dB,  indicating  a  very  weak 
clutter-edge  effect.  (The  weakness  is  due  to  the  fact  that  the  worst 
clutter  in  the  filtered  data  is  the  same  order  of  magnitude  as  the 
NEI.)  Similarly,  the  threshold  sensor  can  be  made  smaller  than  the 
optimum  size  with  very  little  loss  of  signal  detection;  for  example, 
the  loss  is  only  0.2  dB  at  the  critical  sample  size  in  clutter. 
Consequently,  the  best  1-D  threshold  sensor  size,  under  the  stated 
conditions,  is  the  critical  sample  size  in  clutter  because  the  larger 
optimum  size  risks  a  clutter-edge  effect  but  improves  Bignal 
detection  very  little.  A  2-D  threshold  sensor  probably  would  not  be 
advantageous  in  this  case,  even  though  the  onset  of  the  clutter-edge 
effect  is  slightly  below  the  critical  sample  size  in  clear  sky;  the 
reason  is  that  the  weakness  of  the  clutter-edge  effect  allows  the 
critical  sample  size  in  clear  sky  to  be  reached  with  no  significant 
loss  of  signal  detection. 

Finally,  signals  are  almost  always  detected  at  their  peaks, 
although  rarely  a  signal  missed  at  the  peak  is  detected  elsewhere  on 
the  waveform.  Off-peak  detection  is  more  likely  at  lower  detection 
frequencies  but  is  never  very  important.  Specifically,  at  detection 


frequencies  >  0.5,  a  signal  missed  at  the  peak  has  less  than  a  1  in 
80  chance  of  being  detected  elsewhere  on  the  waveform.  Although  this 
result  was  obtained  from  computations  with  LMS-filtered  noise,  it  is 
probably  true  for  any  noise  with  zero  local  mean  and  correlation 
length  about  the  same  as  the  signal  length. 


Appendix.  Summary  Descriptions  of  Data,  Algorithms,  and 
Method  of  Simulation 

(This  appendix,  excepting  Fig.  9,  summarizes  the  detailed 
descriptions  of  the  same  topics  in  the  previous  article.*  The 
summary  is  provided  at  the  request  of  the  editor  and  reviewers  to 
make  this  article  reasonably  self-contained.) 

The  background  data  were  collected,  recorded,  and  digitized  by 

14 

the  General  Electric  Company.  The  optical  sensor  consisted  of  a 
scanning  telescope  with  a  4.0-4. 8  pm  spectral  filter  and  an  InSb 
detector  of  six  photosensitive  elements  in  the  telescope's  focal 
plane.  Each  element  was  0.15  mR  in  azimuth  (the  scan  direction)  and 
0.36  mR  in  elevation;  the  element  width  was  matched  to  the  diameter 
of  the  point-source  image.  Data  collection  took  place  at  43°N,  76°W 
on  22  February  1974  at  5  p.m.  with  the  telescope  pointed  in  a 
westerly  direction  at  an  elevation  of  12-13°.  The  temperature  was 
60°F,  relative  humidity  70%,  and  visibility  25  km  or  more.  The  scene 
contained  clear  sky  (38%  of  the  total  field  of  view)  and  midlevel 
cumulus  clouds  (62%  of  the  field!  with  extremely  bright  edges  backlit 

by  the  setting  sun.  The  range  to  the  clouds  was  10-20  km,  as 

estimated  from  their  type  and  the  conditions  of  observation. 

The  steps  taken  in  recording  and  digitizing  the  detector  outputs 

are  illustrated  in  the  top  half  of  Fig.  7.  Owing  to  the  brightness 
of  the  cloud  edges,  wideband  filters  were  needed  to  limit  the  outputs 
to  the  dynamic  range  (42  dB)  of  the  analog  tape  recorder.  The 
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passband  was  sufficiently  wide  so  that  its  effect  on  the  noise  and 

signal  (the  point-source  output  of  a  detector  element.  Fig.  8)  waa 

negligible  compared  to  the  effect  of  the  LMS  filter  simulated  in  this 

work.  On  the  other  hand,  the  low-frequency  time  constant,  (1/40) 

sec,  was  too  small  relative  to  the  scan  time,  (3/4)  sec,  to  allow 

radiometric  calibration  of  the  data.  The  noise  equivalent  irradiance 

-13  2 

(NEI),  referred  to  the  wideband  output,  was  1.5  x  10  W/cm  . 
Recording,  sampling,  and  digitization  did  not  increase  this  value, 
l.e.,  did  not  add  more  noise  to  the  data.  The  sampling  rate  waa 
about  3.5  per  dwell,  and  the  sampling  time  was  much  less  than  the 
dwell  time  (0.344  msec).  Digital  values  consisted  of  a  sign  bit  and 
eleven  numeric  bits. 

Wideband  output  from  a  typical  cloud  edge  is  plotted  in  Fig.  2A. 
The  maximum  clutter-to-NEI  ratio  is  about  36  dB.  Inspection  of  the 
plot  and  physical  considerations  indicate  that  the  clutter  is  the  sum 
of  two  components:  a  deterministic,  pure  clutter  component  from 
spatial  changes  of  cloud  radiance,  and  a  random  component  from 
photon-arrival  fluctuations  and  charge  fluctuations  in  the  detector 
and  electronics.  (When  the  random  component  exists  alone,  it  is 
often  called  NEI  noise  or  detector  noise.)  Fig.  2B  shows  the 
cloud-edge  clutter  of  Fig.  2A  after  passage  through  the  LMS  filter 
described  below.  The  maximum  clutter-to-NEI  ratio  is  reduced  about 
20  dB,  and  the  clutter  below  158  mR  and  above  165  mR  is  visually 
Indistinguishable  from  LMS-filtered  uniform-sky  noise.  Obviously  the 
LMS  filter  suppresses  the  pure  clutter  component  very  effectively. 


However,  clutter  suppression  alone  is  not  an  adequate  measure  of 
filter  performance;  the  signal-to-clutter  and  signal-to-noise  ratios 
given  below  are  needed  also. 

The  LMS  filter  is  a  linear,  time-invariant,  electronic  filter; 

with  input  from  a  scanning  optical'  sensor,  it  acts  as  a  spatial 

filter.  Its  impulse  response  is  derived  as  follows:  The  point- 

source  signal  and  additive  noise  are  modeled,  respectively,  by  a 

raised  cosine  function  s(t)  and  By  the  first  two  terms  of  a  Maclaurin 

series  n(t)  *  b^  +  b^t  +  -  -  -;  thus  the  sum  of  noise  and  a 

signal  with  arbitrary  amplitude  A  is  m(t)  ■  As(t)  +  (b^  +  b^t) . 

This  model  is  fit  by  least  squares  to  the  filter  output  v(t)  over  one 

period  of  the  raised  cosine,  resulting  in  integral  expressions  for 

the  best-fit  values  of  A,  bg,  and  b^.  The  equation  for  the  best-fit 

A  is  a  convolution  integral  which  yields  the  LMS  filter’s  impulse 
13 

response.  Fig*  8  shows  this  Impulse  response  and  output  from  a 
point-source  input;  Fig.  9  exhibits  the  filter's  amplitude  response 
together  with  the  square  roots  of  the  clutter  power  density  and 
signal  energy  density.  (The  power  density  of  the  uniform-sky  noise 
is  approximately  constant,  to  the  wideband  cut-off  at  5600  Hz.)  The 
.impulse  response  is  shaped  like  the  signal's  mirror  image  but  lowered 
so  that  its  average  value  is  zero.  The  shaping  matches  the  filter  to 
the  signal;  the  symmetry  and  lowering  cause  the  filter  to  eliminate 
the  dc  and  entire  odd  part  of  any  input.  Together  the  shaping  and 
lowering  of  the  impulse  response  result  in  the  LMS  filter  having  the 
orjlmum  ratio  of  peak  signal  to  rms  noise  for  a  point-source  signal 


in  noise  whose  power  spectrum  varies  as  the  reciprocal  of  the  squared 
electrical  frequency.  Such  a  function  models  cloud  clutter  noise. 

The  LMS  filter  is,  therefore,  not  optimum  in  clear  sky,  but  it 
provides  a  good  compromise  between  performance  in  clear  sky  and  in 
cloud  clutter.  For  example,  the  LMS  filter  has  a  maximum  signal— to- 
clutter  ratio  which  is  22  dB  greater  than  that  of  the  low-pass  filter 
with  the  best  possible  signal-to-noise  ratio  in  clear  sky.  The  price 
paid  for  this  large  advantage  in* clutter  is  only  1.3  dB  of  clear-sky 
signal-to-noise  ratio,  relative  to  the  same  low-pass  filter.* 

The  adaptive  threshold  sensor  of  Fig.  1  is  discussed  in  section 
II  where  it  is  explained  that  this  sensor  makes  a  statistical  test 
-for  signals  if  the  noise  mean  of  the  filter  output  is  zero  locally 
and  globally  (over  extents  £  the  sensor's  length).  Uniform-sky  noise 
has  this  type  of  mean,  as  does  the  cloud-clutter  output  of  a  filter 
with  a  high  cut-on  frequency  like  that  of  the  LMS  filter  (Fig.  2B) . 
With  a  low  cut-on  frequency  like  that  of  the  wideband  filter,  the 
clutter  mean  is  zero  globally  but  not  locally  (Fig.  2A).  With  a 
low-pass  filter,  whose  cut-on  frequency' is  zero,  the  clutter  mean  is 
non-zero  both  locally  and.  globally.  The  sensor  of  Fig.  1  cannot 
.operate  statistically  with  the  clutter  output  of  a  filter  having  a 
low  or  zero  cut-on  frequency,  because  it  cannot  estimate  local  values 
of  the  mean  p  and  mean  absolute  deviation  i)  for  setting  the  threshold 
(p  *  Kn).  Rather,  the  sensor  estimates  the  local  noise  magnitude, 
which  is  large,  and  sets  the  threshold  at  a  level  several  times 
greater,  virtually  turning  off  the  surveillance  set  in  clutter. 


A  filter  should  condition  the  noise  so  that  the  threshold  sensor  can 
function  statistically.  The  LMS  filter  does  this  by  reducing  the 
correlation  length  in  clutter  so  that  the  noise  average  is  zero 
within  the  sensor  frame. 

The  simulation  was  carried  out  with  a  digital  computer  using 
programmed  algorithms  for  the  LMS  filter  and  adaptive  threshold 
sensor.  Fig.  1  is  effectively  a  flow  chart  from  which  the  reader  can 
program  the  sensor  algorithm,  given  also  that  the  shutter  prevents 
the  wings  of  a  point-source  signal  centered  on  the  signal  window  from 
distorting  the  estimate  of  the  local  noise  dispersion.  An 
approximation  described  in  section  11(C)  can  be  used  to  simplify  the 
sensor  algorithm  with  negligible  error;  The  principles  underlying 
any  digital  filter  algorithm  in  both- the- time  and  frequency  domains 
are  discussed  in  many  texts,  edited  books,  and  review  articles; 
reference  15  gives  an  especially  lucid,  introductory  account; 
references  1  and  13  describe  the  time-domain  LMS  algorithm  used  in 
this  work.  After  studying  reference  15,  a  reader  with  no  prior 
knowledge  of  digital  filter  principles  can  program  the  LMS  algorithm 
with  understanding.  Regrettably,  the  principles  cannot  be  briefly 
.and  intelligibly  summarized  for  him.  However,  if  he  is  content  to 
program  the  algorithm  without  understanding,  two  facts  are  all  he 
needs:  Each  LMS  filter  output  is  a  weighted  sum  of  the  present  and 
the  M  consecutive  past  inputs,  and  the  (M+l)  normalized  weights  can 
be  obtained  from  the  recipes  given  in  reference  1,  section  IV. 

The  main  steps  of  the  simulation  are  indicated  in  the  lower 
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half  of  Fig.  7.  After  LMS-filtering  of  the  background  data,  the 
number  of  false  alarms  was  determined  as  a  function  of  the  threshold 
parameter  K.  This  was  done  by  calculating  the  ratios  K  “  xw/<|x^}> 
and  ordering  them  in  statistical  frequency  tables.  The  value  of  K 
which  allowed  one  false  alarm  in  All  the  data  was  then  read  from  the 
tables  and  given  to  the  threshold  sensor  algorithm.  Aa  each  noise 
value  entered  the  signal  window  it  was  added  to  a  simulated  signal 
amplitude,  and  the  sum  was  tested  against  the  threshold  K<|x^|>. 
(Because  the  filter  is  linear,  simulated  signals  can  be  added  to  the 
noise  either  before  or  after  the  filter,  provided  the  signal  waveform 
appropriate  for  each  case  is  used.)'  The  simulated  signal  amplitude 
was  then  increased  and  the  process  repeated  to  determine  the  curves 

3 

of  detection  frequency  in  Fig.  5.  This  procedure  assumes  that  a 
signal  is  detected  only  at  its  peak,  an  approximation  which  introduces 
<  2.5%  error  in  detection  frequencies  t  0.5. 
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Figure  Captions 

Figure  1.  Schematic  diagram  of  an  adaptive  threshold  sensor.  The 

parameter  K  is  fixed,  but  the  threshold  T  varies  with  the 

estimate  <|x^|>  of  the' local  noise  mean-absolute- 

deviation.  The  value  of  K  required  to  prevent  a  false 

alarm  from  window  value  x  is  x  /<|x_|>. 

w  w  f 

Figure  2.  Plots  of  backlit  cloud-edge  clutter  at  the  outputs  of  the 
(A)  wide  band  and  (B)  LMS  filters.  The  wide  band  output 
provided  the  noise  data  used  in  this  work.  The  NEI  unit 
of  the  ordinate  is  the  standard  deviation  (and  rros 
voltage)  of  the  clear-sky  noise  from  the  wide  band  filter. 

Figure  3.  .Sampling  standard  deviations  (s.d.s)  of  sensor  frame- 

averages  n  ■  <|x^|>  as  functions  of  sample  (frame)  size, 
noise  type,  and  segment  length.  The  noise  means  p  are 
zero  locally  and  globally,  so  that  noise  magnitude  }x| 
and  absolute  deviation  |x  -  p|  are  identical,  and  their 
means  are  estimated  by  f}  locally  and  globally.  The 
sample  sizes  are  8,  18,  36,  80,  144,  240,  340,  440,  540, 
and  640  data.  Noise  data  of  each  type  are  divided  into 
segments  of  the  specified  lengths.  There  are  nine 
segments  in  each  of  sets  A-D,  and  ten  in  each  of  sets  E-F. 
Variances  are  calculated  about  q  ■  <|x|>,  the  average 
magnitude  and  average  absolute  deviation  of  all  data  in  a 
set.  Bars  show  the  ranges  of  the  sampling  s.d.s  of  the 
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Figure  4 


Figure  Captions  (cont.) 

segments;  dots  indicate  the  sampling  s.d.  of  the 
population  (all  segments  of  the  set).  A  population 
sampling  variance  is  equal  to  the  average  of  the 
associated  segment  sampling  variances.  Sampling  s.d .s 
0-  are  divided  by  the  associated  population  s.d.  0  to 
facilitate  comparison  of  results  from  different  data  sets. 
The  diagonals  represent  the  logarithmic  form  of 
equation  (1),  vhich  is  the  relation  between  sample  size 
and  sampling  s.d.  given  by  the  theorem  on  the  sampling 
distribution  of  the  mean.  The  lines  through  the  dots  of 
data  sets  C-F  are  drawn  for  best  visual  fits;  departures 
from  the  diagonals  are  due  to  non-random  sampling, 
correlation,  and  non-stationarity  of  the  noise. 

Sampling  standard  deviation  (s.d.)  and  threshold  parameter 
K  as  functions  of  threshold-sensor-sample  size  and  noise 
type.  The  sample  sizes  are' 36,  80,  144,  240,  and  340 
data.  (A)  Population  sampling  s.d.s  for  the  data  sets  of 
Figs.  IE  and  IF',  each  set  contains  10  segments. 

(B)  Threshold  parameters  for  one  false  alarm  (f.a.) 
determined  from  the  two  data  set3  of  reference  1,  each  of 
which  contains  114  segments. 


Figure  5.  Frequencies  of  detecting  simulated  signals  in  LMS-filtered 
natural  background  noise  with  threshold  sensors  of 


Figure  Captions  (cont.) 


Figure  6 


Figure  7 


Figure  8 


different  sizes.  These  results  are  calculated  from  the 
two  data  sets  used  in  reference  1;  each  set  contains  114 
segments.  The  difference  between  signals  and  pseudo 
signals  is  explained  in  section  11(D).  Results  for  a 
frame  size  of  144  data  are  in  reference  1;  they  are 
hardly  distinguishable  from  the  results  for  a  frame  size 
of  340  data.  The  threshold  parameters  for  one  false  alarm 
(f.a.)  are  plotted  in  the  top  curve  of  Fig.  4B. 

Equally  detected  signal  strength  as  a  function  of 
threshold-sensor  size,  noise  type,  and  detection 
.  frequency.  Simulated  signal  strengths  for  these  plots 
were  read  from  the  curves  in  Fig.  5.  The  curves  here  show 
simulated  signal  strengths  required  to  achieve  the 
indicated  detection  frequencies  with  threshold  sensors  of 
different  sizes.  The  difference  between  signals  and 
pseudo  signals  is  explained 'in  section  11(D). 

Steps  taken  in  collecting  the  background  data  and 
performing  the  biraulation. 

Time-domain  properties  of  the  point-source  signal  (input) 
and  the  LMS  filter.  The  output  is  the  convolution  of  the 
input  and  the  LMS  impulse  response.  The  tick  marks  are 
one  dwell  time  apart. 
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FACTORS  AFFECTING  BB  PULSE  SHAPE 


1.  Orientation  of  BB  relative  to  scan 

2.  Vertical  positioning  of  BB  image  on  detector  array 

3.  Flux  distribution  in  defocused  image 
A.  Alinement  of  optics  and  detectors 

5,  Response  of  electronics 

6.  BB  temperature  distribution 


FIELD  CALIBRATION  ALTERNATIVES 


1.  Distant  extended  source 

2.  Close  source 

a.  Extended 

b.  Point  (Jones  method) 

3.  "Collimator"  (spectacles) 
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EXTENDED  SOURCE  WISH  LIST 


i 


Three  panels  with  temperature  ranges 

a.  190  to  240  K  (-80  to  -30  C) 

b.  240  to  290  K  (-30  to  20  C) 

c.  290  to  340  K  (20  to  60  C) 

Panel  size:  18  x  18  inches 


COMMERCIALLY  AVAILABLE  SOURCE 

Barnes  11-140T 
Single  panel,  12  x  12  inches 
Temperature  range:  ambient  to  500  K  (230  C) 


2 


Introduction 

The  optical  sensors  used  in  the. Background  Measurements  and  Ana¬ 
lysis  Program  (BMAP)  are  to  be  calibrated  in  the  laboratory  so  that  data 
collected  in  the  field  can  be  assigned  absolute  radiance  values.  Internal 
references  in  the  sensors  can  be  used  to  link  their  responses  in  the  field  to 
the  laboratory  calibrations,  but  this  provides  only  one  calibration  point,  and 
at  least  two,  preferably  three,  are  needed  to  estimate  both  reference  radiance 
and  responsivity.  To  strengthen  the  link  between  laboratory  calibration  and 
field  response,  a  second  calibration  in  the  field  is  desirable. 

This  note  considers  the  questions,  "What  is  the  purpose  of  a  field  cali¬ 
bration  (already  answered),  how  should  it  be  performed,  and  what  information 
must  it  provide?"  Concomitantly,  the  first  steps  toward  reducing  the  field 
calibrations  made  during  the  background  measurements  at  Montauk  Point  are 
described. 

To  begin,  the  possible  methods  of  calibration  are  reviewed.  "Calibra¬ 
tions  are  categorized  into  two  groups:  those  in  which  the  radiometer  field 
stop  is  partially  illuminated  by  the  calibration  source;  and  those  in  which 
the  stop  is  totally  filled  with  uniform  calibration  radiation."  (The  Infra¬ 
red  Handbook,  p.  20-13.)  The  first  type  of  calibration  is  impractical  for 
field  work  with  a  multichannel  detector  because  the  source  image  must  be 
carefully  positioned  several  times  to  calibrate  all  channels.  BMAP  field 
calibrations  must  be  of  the  second  type.  "The  latter  can  be  realized  with  a 
so-called  extended  source,  either  distant  --  or  close  --,  by  a  collimator  — 
or  by  a  proper  placement  of  a  point  source  very  close  to  the  entrance  aper¬ 
ture  (the  so-called  'Jones  method'  --)."  (ibid.) 


The  Jones  method  is  attractive  in  several  ways:  it  can  be  used  both  in 
flight  and  on  the  ground;  the  source. is  easily  accessible  to  the  operators; 
it  does  not  have  to  be  carefully  positioned;  and  atmospheric  absorption  is 
negligible  so  that  the  spectral  distribution  can  be  the  same  for  the  field 
and  laboratory  calibration.  Two  disadvantages  are  the  need  for  very  small 
standard  sources  and  means  of  positioning  them  so  that  the  mounts  do  not 
radiate  significantly  into  the  sensor.  Also,  because  the  sources  are  not 
imaged  on  the  detectors,  they  cannot  be  exposed  simultaneously  to  the  sensor. 
Furthermore,  radiation  received  by  a  detector  element  passes  through  only  a 
smallpart  of  the  lens,  not  all  of  it  as  in  the  field  measurements  and 
laboratory  calibrations.  The  same  comments  apply  to  calibration  with  a  close 
extended  source  except  the  remarks  concerning  source  size  and  aperture 
differences.  Chiefly  because  two  sources  cannot  be  used  simultaneously, 
close-source  calibrations  are  not  completely  satisfactory,  but  neither  are 
the  alternatives,  so  the  close  source  methods  cannot  be  ruled  out. 

Use  of  a  collimator  has  the  disadvantage  that  two  would  be  required  to 
deploy  two  sources.  This  leaves  the  distant  extended  source  method.  Since 
it  was  used  at  Montauk  Point,  it  will  be  explored  more  fully  with  reference 
to  conditions  there.  These  conditions  result  in  a  defocused  image  on  the 
detector.  .A  problem  concerning  the  size  of  this  image  is  stated,  a  solution 
is  given,  and  the  implications  of  the  solution  are  discussed.  Finally  the 
shapes  of  the  calibration  outputs  recorded  on  Montauk  Point  are  considered, 
and  a  procedure  is  described  for  analyzing  them  to  obtain  calibrations. 

Problem 

The  radiometer  is  focused  at  infinity  and  then  pointed  at  a  blackbody 
(BB)  placed  a  large  but  finite  distance  from  the  lens.  Thus,  the  blackbody 
image  lies  behind  the  detector  surface,  which  is  in  the  focal  plane.  There 


are  two  optical  systems,  long  wave  (LW)  and  mid  wave  (MW),  with  parallel  axes 
4(5/8)  inches  apart  in  the  vertical  plane.  How  does  the  size  of  the  defocused 


"image"  on  the  detector  compare  to  the  size  of  the  true  (focused)  image? 
Assumptions 

The  solution  is  based  on  the  thin,  lens  equations  of  geometric  optics. 
Aberrations  are  ignored.  The  center  of  the  88  is  assumed  to  be  on  an  optic 
axis.  (When  the  B8  is  far  away  its  center  is  more  likely  to  be  on  or  near 
the  mid-rline  between  the  optic  axes,  but  this  results  in  negligible  error, 
as  verified  by  computation.)  The  true  image  is  required  to  fill  or  overfill 
the  instantaneous  field  of  view  (IFOV)  in  the  scan  (horizontal)  direction  and 
the  total  field  of  view  (TFOV)  in  the  cross-scan  (vertical)  direction. 

Notation  and  Optical  Parameters  at  Montauk  Point 

f  *  focal  length  of  lens,  6  inches 

r.  »  lens  (or  aperture)  radius,  2  inches 

Horizontal  IFOV  =  (1/3)  mR 

Vert.ical  TFOV  =  5(1/3)  mR 

g  =  object  (B8)  distance,  840  inches 

G  =  object  (BB)  size,  12  x  12  inches. 

i  =  true  (focused)  image  distance 

I  =  -true  image  size 

Subscripts  h  and  v  denote  the  horizontal  and  vertical  directions 
A  denotes  a  change  in  any  quantity. 

Figure  1  shows  the  optical  geometry. 


Solution 


M 

I 


From  Figure  1  and  the  thin  lens  equations: 


.  2  A i  tan  9  =  2(i-f)  fr-(I/2) 


I 


■  nr-) 


------  or  i  =  gf/(g-f) 

i  f  .  g 


^h,v  /Gh,v  “  i/9  or  lh,v/f  "  Gh,v/(9"0 

Eliminating  image  parameters  from  the  right  hand  side  of  the  equation 
for  (a I /I )  gives: 

/AiV  ,  .  !  r  2,r(fl-f.L  .  {\ 

\  I  /h,y  9  LfGh.v  j 


_  Zr 

(l  ‘  ^ 

+  f  =  2r 

f  f  Zr 

--  ll 

Gh,v 

9  Ghs  v 

9  \Gh,v 

) 

Since  the  true  image  must  fill  or  overfill  the  field  of  view,  the  frac¬ 
tional  enlargement  (Al/IJ^v  the  defocused  image  is  constrained  by  the 
two  relations, 

Ih,v/f  «  Gh,v/(g-f)  l  IF0Vh  or  TF0Vv. 

Discussion  of  BB  Image  Size 

This  solution  can  be  checked  against  measurements  made  at  Montauk  Point 
From  the  optical  parameters  given  above: 

( A I / 1 )h>  v  =  (1/3)  +  (2f/3g)  =  0.3381 


i  =  6.04317  in.  Ih>v  =  0.08633  in. 

Ih, V/F  =  14.39  rrR  >  IF0Vh  or  TF0Vv 
Alh,v  =  4*87  nR 

Defocused  image  size  =  (I  A I ) v  =  19.26  mR 


If  this  defocused  image  were  scanned  across  a  detector  element  (1/3)  mR  wide, 
an  output  pulse  with  a  base  19.6  mR  wide  would  be  produced.  Figures  (3)  -  (6) 
show  pulses  recorded  at  Montauk  Point;  their  bases  are  about  19.2  mR  wide,  close 
enough  to  the  calculated  width  to  validate  the  solution.  For  comparison. 

Figure  2  shows  an  idealized  pulse  that  would  be  produced  by  the  true  image. 

(Further  explanation  of  Figures  (3)  -  (6)  is  in  order.  The  first  112 
samples  of  the  plots  display  the  internal  reference  pulse  and  the  clamped 
amplifier  output  with  the  detector  viewing  itself.  There  are  two  intensity 
scales.  The  one  on  the  left  is  for  the  part  of  the  figure  leftward  of  the 
tallest  vertical  line  in  Figure  3;  the  right-hand  scale  is  for  the  rest  of  the 
figure.  The  saturation  intensity  is  4096,  so  the  internal  reference  source 
almost  saturates  the  IW  output.  The  dashed  spikes  are  recording  errors;  the 
second  horizontal  axis  and  ticks  are  related  to  these  errors.  The  errors  are 
of  no  interest  here  because  they  will  be  corrected  before  the  calibration  is 
derived. ) 

It  is  useful  to  consider  how  (aI/1)  depends  on  the  optical  parameters. 

In  eq.  (1)  ordinarily  l»(f/g)«(2r/G);  consequently  the  fractional  enlarge¬ 
ment  of  the  defocused  image  is  determined  mainly  by  the  ratio  (2r/G).  With 
a  finite  object  both  A I  and  I  approach  zero  as  g  approaches  infinity,  so 
(AI/-I)  is'  indeterminate.  However,  using  VHospital's  rule  one  finds 

lim  (AI/I  )g  -•  =  (2r/G). 

This  agrees  with  the  approximation,  (Al/I)  =s  (2r/G),  for  large  g.  Differen¬ 
tiating  eq.  (1)  with  respect  to  object  distance  gives 
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FIGURE  3:  Pulse  from  long-wave  channel  2  (low  elevation) 
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FIGURE  2:  Ideal  pulse  produced  by  true  image 


Schematic  optical  geometry 


and,  if  so,  whether  corrections  can  be  made  which  allow  true  pulse  heights  to 
be  computed. 

If  true  pulse  heights  can  be  determined,  the  next  step  is  to  calculate 
the  external  BB  radiance  at  the  lens  using  the  LOWTRAN  atmospheric  transmission 
and  emission  model.  This  radiance  and- the  true  pulse  heights  constitute  one 
calibration  point.  The  internal  BB  radiance  and  the  heights  of  the  internal 
reference  pulses  provide  another,  but  not  directly  because  the  internal  BB 
radiance  is  relayed  to  the  detector  by  one  optical  system,  the  external  BB 
radiance  and  the  scene  radiance  by  another,  and  the  two  optical  systems  trans¬ 
mit  and  emit  differently.  Thus,  the  internal  BB  radiance  must  be  corrected 
once  to  obtain  radiance  on  the  detector,  and  a  second  time  to  obtain  the 
equivalent  radiance  at  the  sensor  aperture.  The  sensor  response  to  scene 
radiance  can  then  be  derived  from  the  calibration  points. 

Note  that  this  procedure  gives  two  values  for  responsivity  but  none  for 
reference  radiation.  Perhaps  the  latter  can  be  estimated  from  the  closed 
cover. scans;  this  remains  to  be  seen. 


and  centers  of  the  pulses  had  different  but  regularly  varying  widths.  Unfortu¬ 
nately,  the  pulses  on  some  channels  at  other  elevations  look  very  different, 
e.g.,  compare  Figure  3  with  Figure  4,  and  Figure  5  with  Figure  6.  If  the 
amplifier  for  MW  channel  1  was  working  properly,  the  pulse  shape  in  Figure  5 
must  have  resulted  from  some  combination  of  non-uniform  flux  distribution  and 
maladjustment  of  factors  (2)  r  (4)  above.  At  this  time  only  one  set  of  LW 
calibration  pulses  and  two  of  MW  pulses  have  been  plotted.  This  is  not  enough 
to  determine  more  precisely  the  cause  of  the  deformed  pulses  in  Figures  (4) 
and  (5). 

For  calibration  one  needs  to  know  the  heights  of  the  pulses  that  would  be 
produced  by  the  true  image  crossing  the  detector  elements.  Assuming  conditions 

(1)  -  (3)  and  (6)  are  ideal,  the  true  heights  can  be  computed  from 

* 

(a)  the  heights  and  areas  of  the  actual  pulses, 

(b)  the  pixel  size, 

(c)  the  optical  parameters, 

(d)  the  vertical  position  of  the  defocused  image  on  the  detector  array,  and 

(e)  the  flux  distribution  in  the  defocused  image. 

The  flux  distribution  can  be  calculated  if  the  radiance  is  uniform  over  the  true 
image  area  and  isotropic  over  the  range  of  solid  angle  within  the  illuminated 
volume  of  the  optical  system.  In  that  case  the  calculation  amounts  to  computing 
the  projected  true-image  area  which  receives  radiation  from  a  point  of  the 
defocused  image.  If  other  conditions  are  ideal  as  assumed,  the  vertical 
position  of  the  defocused  image  can  be  inferred  from  the  flux  distribution  and 
the  actual  pulse  "shapes.  At  Montauk  Point,  however,  one  or  more  of  the  other 
conditions  may  not  have  been  ideal.  Plots  of  all  calibration  scans  at  Montauk 
must  be  Inspected  to  determine  whether  any  conditions  were  in  fact  non-ideal 


(5)  the  flux  distribution  in  the  defocused  image,  and 

(6)  the  response  of  the  electronics. 

The  first  factor  can  be  ignored  for  present  purposes;  the  blackbody 
temperature  is  known  to  be  uniform  in  comparison  to  the  effects  of  other 
factors.  The  amplifiers  for  the  LV1  array  are  ac  coupled  with  a  0.5  Hz,  3  dB 
cut-on  frequency;  those  for  the  MW  array  are  dc  coupled.  The  3  dB  cut-off 
frequency  for  both  sets  of  amplifiers  is  1  kHz  (1  ms  time  constant).  The 
spinning  mirrors  for  the  internal  reference  source  are  shaped  so  that  full 
response  to  its  radiation  is  achieved,  but  only  momentarily.  Before  each 
scan  the  amplifier  outputs  are  clamped  with  the  detectors  viewing  themselves. 

Assume  conditions  (1)  -  (3)  above  are  ideal  and  the  true  image  is 
focused  on  and  vertically  centered  on  the  detector  array.  Then  in  principle 
the  B8  pulse  from  the  dc  coupled  amplifiers  appears  as  in  Figure  2.  The  pulse 
from  the  ac  coupled  amplifiers  is  the  sane  except  that  its  top  droops  slightly 
and  there  is  a  slight  undershoot  after  the  pulse.  Some  of  the  actual  pulses 
(Figures  3  and  6)  resemble  the  ideal  pulses,  but  they  have  wider  bases,  wider 
shoulders,  and  narrower  tops.  The  wider  bases  result  from  the  increased  width 
of  the  defocused  image  as  explained..  The  wider  shoulders  and  narrower  tops 
can  be  explained  by  a  non-uniform  flux  distribution  across  the  width  of  the 
defocused  image.  The  droop  of  the  pulse  top  in  Figure  3  is  only  4  per  cent 
greater  than  that  expected  from  the  cut-on  frequency  of  the  ac  coupled  am¬ 
plifiers,  but  the  droop  in  Figure  4  is  more  than  twice  as  great,  and  the  top 
has  an  unexpected  upward  bow. 

If  the  puVses  on  all  channels  had  the  same  base  widths  and  appearance  as 
those  in  Figures  3  and  6,  the  results  could  be  explained  by  the  size  and 
non-uniform  flux  distribution  of  the  defocused  image,  even  if  the  shoulders 


at  which  distance  (Al/I)  =  0.3369.  At  the  actual  distance  of  70  ft.  (Al/I)  = 

0.3381,  so  at  Montauk  Point  the  BB  was  about  as  far  from  the  radiometer  as  it 
profitably  could  have  been  from  the  standpoint  of  image  size.  (From  the 
standpoint  of  flux  uniformity  in  the  image,  moving  the  BB  farther  away  might 
or  might  not  have  improved  matters.  The  equations  for  image  size  say  nothing 
about  this;  it  must  be  studied  separately.)  With  respect  to  image  size,  there 
are  only  three  ways  to  improve  upon  the  Montauk  Point  conditions:  use  a 
"collimator,"  use  a  larger  blackbody,  or  reduce  the  sensor  aperture  with  a 
removable  stop  for  calibration.  ,  t 

-  -  -  -  (tiowevcr j  be  possible  to  cu.ro  the  nJionuetef' s 

■far- sijhtedHcss  with  spectacles^  a  related  technique.) 

The  disadvantage  of  a  collimator  has  already  been  noted. ^  With  a  larger  BB, 
temperature  uniformity  is  harder  to  achieve.  In  addition  a  larger  BB  may  not 
be  commercially  available  and  is  probably  expensive  if  it  is  available.  Cost 
is  all  the  more  critical  because  at  least  two  sources  are  needed.  If  the 
aperture  is  reduced  for  field  calibrations,  it  will  be  different  from  that  for 
laboratory  calibrations  and  field  measurements.  (Recall  the  similar  disad¬ 
vantage  of  the  Jones  method.)  On  the  other  hand,  the  purpose  of  the  field 
calibration  is  to  link  the  radiometer's  field  response  to  the  laboratory 
calibration.  An  aperture  stop  should  be  used  if,  by  reducing  defocusing 
effects,  it  strengthens  the  link.  This  would  have  to  be  determined  experi¬ 
mentally.  The  stop  could  be  a  small  hole  in  the  aperture  cover. 

Discussion  of  BB  Pulse  Shape  and  Calibration  Analysis 
The  shapes  of  the  BB  pulses  are  affected  by 

(1)  the  BB  temperature  distribution, 

(2)  orientation  of  the  BB  relative  to  the  scan  direction, 

(3)  alignment  of  the  optics  and  detectors, 

(4)  vertical  positioning  of  the  BB  image  on  the  detector  array. 


This  is  positive  for  G  <  2r,  zero  for  G  =  2r,  and  negative  for  G  >  2R.  Thus, 
as  the  BB  is  moved  away  from  the  radiometer,  the  ratio  (aI/I) 


increases  to  (2r/G)  if  G  <  2r. 
remains  equal  to  (2r/G)  if  G  *  2r. 
decreases  to  (2r/G)  if  G  >  2r. 


Obviously  there  are  conditions  for  which  (aI/I)  =  0.  In  this  case 

v\fM)- 

If  G  <  2r  the  object  must  be  inside  the  focus  and  the  image  virtual  for 
(aI/I)  -  0.  This  result  is  related  to  both  the  close  extended  source  and 
Jones  calibration  methods.  Its  implications  will  not  be  considered  here. 

If  G  >  2r  a  virtual  object  is  required  for  (aI/I)  =  0.  This  is  achieved  with 
a  collimator.  (Here  "collimator"  is  a  misnomer  since  the  device  forms  an 
image,  not  a  parallel  beam.) 

With  a  distant  extended  source  G  >  2r,~sd  moving  the  BB  away  from  the 
sensor  reduces  the  fractional  enlargement  of  the  defocused  image,  but  very 
slowly  at  large  object  distances  and  not  below  (2r/G).  In  any  case  the 
object  distance  is  limited  by  the  object  size  since  the  true  image  must 
fill  or  overfill  the  field  of  view,  perferably  by  two  or  more  times  to  make 
positioning  of  the  image  on  the  detector  less  critical.  In  order  to  satisfy 
this  overfill  requirement,  the  distance  to  the  BB  at  Montauk  Point  could  not 
have  exceeded 


g  =  Gv 
2  ( T  F  0  V  v ) 


•  I 


+  f  =  1311  in.  =  94  ft 


Sr/ia  II 


STANDARD  SOURCE 


SPECIFICATIONS 

Area:  10  in.  horizontal  (minimum) 

10  in.  vertical  (minimum) 

Temperature  Range:  -20  to  50  °C  (253-323  °K) 

MODIFIABLE  TO  -50  °C  (223  °K) 

Uniformity:  0.1  °C  desired  (some  non-uniformity 
near  edges  allowed) 

Temperature  Readout  Resolution:  0.1  °C 
Analog  Controller  Acceptable 


MAKERS 


Advanced  Kinetics 
Barnes  Engineering 
Block  Engineering 


DBA  Systems 
Electro-Optical  Inds. 
Eppley  Laboratory 
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TYPES  OF  CALIBRATIONS 


NEI  and  Spatial  Response 
(Point  and  Bar  Source  Scans) 


Radiometric  Response 
(Near  Extended  Source  Scans) 


LABORATORY  CALIBRATION  5ET-U.P 


Blackbody  sou-rce 


LABORATORY  CALIBRATION  PROCEDURE 


1.  Point  Source  Scans 

a.  Center  point-source  image  vertically  on  pixel. 

b.  Record  4  frames, 

C.  Do  A-B  FOR  ALL  CHANNELS  OF  DETECTOR  ARRAY. 

D.  Do  A-C  FOR  BOTH  ARRAYS. 

2.  Bar  Source  Scans 

a.  Center  bar  target  vertically  on  detector  array. 

b.  Record  4  frames  in  alternating  mode  (2  frames 

FROM  EACH  DETECTOR)  , 

C.  Do  A-B  WITH  BARS  AT  10,  20,  30,  40,  50  °C. 

3.  Near  Extended  Source  Scans 

a.  Remove  bar  apertures  from  source. 

b.  Place  source  close- to  telescope  objective  lens, 

c.  Record  4  frames, 

.  D.  Do  C  WITH  SOURCE  AT  10,  20,  30,  40,  50  °C. 


4.  Do  1-4  WITH  RADIOMETER  AT  20,  25,  30,  35,  40  °C. 


MW  POINT  SOURCE  RESULTS 
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LW  POINT  SOURCE  RESULTS 


FIELD  CALIBRATION  PROCEDURE 


1.  Point  Source  Scans 

a.  Center  image  from  portable  collimator 

VERTICALLY  ON  PIXEL, 

b. .  Record  A  frames, 

c,  Do  a-b  FOR  CHANNELS  2.  8,  15  OF  detector  array. 

d,  Do  a-c  for  eoth  arrays, 


2,  No  BAR  SOURCE  SCANS, 


3,  Near  Extended  Source  Scans 

a.  Same  procedure  as. in  laboratory. 

b.  Record  scans  at  end  of  each  data  tape. 


CALIBRATION  ISSUES  &  PROBLEMS 


POINT  RESPONSE 

New  MW  detector  geometry 

Vertical  registration  of  MW  &  LW  detectors 

Dependence  of  point- image  size  &  shape  on  temperature 

Dependence  of  cross-channel  coupling  on  temperature 

Must  10"  paraboloid  be  centered  on  radiometer  telescope  for 
point-response  measurements? 

Comparisons  between  lab  &  field  measurements 

(aTrl) 


CALIBRATION  ISSUES  &  PROBLEMS 


RADIOMETRIC  RESPONSE 

Calibration  equation:  (v-v0)/(vH-v  )  =  aR  +  b 

LW  HOT  REFERENCE  SATURATES  AT  RADIOMETER  TEMPS  >~30  °C. 

DAS  PROBLEMS:  INTERNAL  HOT  S/OR  COLD  REFERENCES  MISSING  IN 
SOME  CASES 

CAN  RADIOMETRIC  CALIBRATIONS  BE  MADE  IN  THE  FIELD?  (YES) 


(NRL) 


CALIBRATION  ISSUES  &  PROBLEMS 


MISCELLANEOUS 

IS  CHANNEL  12  OF  LW  ARRAY  FAILING? 

Must  radiometer  temperature  be  limited? 

Point-image  distortion  at  high  temperature? 

Increased  cross-channel  coupling  at  high  temperature? 

LW  COLD  CLAMP  MAY  NOT  WORK  AT  HIGH  RADIOMETER  &  CALIBRA¬ 
TION  TEMPERATURES. 

Hot  REFERENCE  PULSE  MAY  NOT  BE  IDENTIFIABLE  AT  HIGH 
RADIOMETER  &  CALIBRATION  TEMPERATURES. 


Data  storage  ("block")  errors 

Data  retrieval  ("bit")  errors 

Slippage  of  data  acquisition  window 
or  internal  reference  mirrors 
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BLOCK  DIAG-RAM :  RAYTHEON  DAS 


Cold  and  hot 
re  -fe,  recces 


POINT  DATA  AUG  23,  1984 

TARGET  TEMPERATURE:  227.1  C  OPTICS  TEMPERATURE:  69.0  F 
3  DAYS  11  HOURS  49  MINUTES  55.8  SECONDS 
CHANNEL  3 

LW  Kolr  /'e f e rence  Lowest 
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MIDWAVE  RADIOMETER  CALIBRATION 


The  radiometer's  response  to  monochromatic  radiation  will  be 
derived  first,  then  the  analysis  will  be  extended  to  broadband 
radiation.  This  will  avoid  hiding'  the  essential  features  of  the 
response  in  the  fog  of  integration  over  the  spectral  band. 

A.  Monochromatic  Response 

Assume  there  is  an  analytic  relation  between  the  detector  output 
voltage  v  and  the  electrical  (E),  optical  (0),  and  thermal  (T) 
parameters  of  the  instrument  and  environment. 

v  *  f(E,0,T)  (1) 

If  all  parameters  except  radiance  r  on  the  detector  are  fixed, 

v  =  g( r) .  (2) 

This  relation  can  be  expanded  in  a  Taylor  series 

v  -  g(0)  +  g ' ( 0) r  +  |g"(0)r2  +  - .  (3) 

The  coefficients  of  the  series  are  functions  of  all  the  system 
parameters  except  r  and,  therefore,  are  constants  only  if  all 
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parameters  except  radiance  on  the  detector  are  fixed.  It  will  be 
assumed  that  the  first  two  terms  of  .the  series  adequately  represent 
the  voltage-radiance  relation 

v  =  g(0)  +  g' (0)r .  (4) 

The  detector  output  is  dc  coupled  to  an  amplifier  which 
multiplies  it  by  a  factor  a  and  adds  a  voltage  B,  so  that  the 
voltage-radiance  relation  for  the  whole  radiometer  becomes 

v  =  [c£g'(0)]r  +  [ag(0)  +  &].  (5 

The  factors  a  and  B  are  constant  only  if  the  electrical  and  thermal 
parameters  of  the  system  are  fixed.  The  added  voltage  is  chosen  so 
that  the  radiometer  output  is  v^  when  the  detector  is  completely 
covered  by  a  nearby  mirror  which  rotates  through  the  field  of  view. 
During  the  time  of  complete  coverage,  the  radiance  r^  on  the  detector 
comes  from  emission  by  the  mirror  and  reflection  of  the  detector's 
own  radiation  back  to  it.  Little  or  no  stray  radiation  from  the 
radiometer's  interior  reaches  the  detector  at  this  time.  The 
voltage-radiance  relation  with  the  cold  reference  mirror  in  place  is 

v  =  (ag'(O)Jr  +  (ag(0)  +  Bj.  (6) 

o  o 

The  cold  reference  mirror  is  followed  through  the  field  of  view  by 


another  mirror  mounted  on  the  same  motor  shaft.  This  mirror  reflects 


radiation  from  an  internal  refetence  source  onto  the  detector.  With 
the  hot  reference  mirror  in  place,  most  of  the  radiance  r^  on  the 
detector  comes  from  the  internal  source  and  the  optics  which  transfer 
its  radiation  to  the  detector.  However,  radiation  from  the 
radiometer’s  interior  is  also  reflected  to  the  detector  at  this  time. 
The  peak  hot  reference  voltage  is 

v^  «  [ag'(0)]rh  +  (etg(O)  +  & ) .  (7) 

The  difference  between  the  hot  and  cold  reference  voltages  yields  a 
relation  for  the  amplification  factor  [ctg'(O)]. 

v,  -  v  «=  (r,  -  r  ){ag’(0)]  (8) 

ho  ho 

The  radiometer  is  calibrated  by  measuring  the  voltage  output 
with  the  total  field  of  view  overfilled  by  an  extended  source  of 
radiance  R  placed  close  to  the  telescope's  objective  lens.  Under 
these  conditions  the  radiance  on  the  detector  is  (tR  +  r) ,  where  r 
accounts  for  losses  of  source  radiation  in  the  main  optics,  and  r  for 
stray  radiation  from  these  optics  and  the  radiometer's  interior. 

With  the  standard  source  in  view,  the  voltage-radiance  relation  is 

(9a) 


v  =  [ ag ' (0) ] ( tR  +  r)  +  [ag(0)  +  8] 

=  [xag'(0)]R  +  (ag'(O)r  +  ag(0)  +  8). 


(9b) 


Equation  (9b)  is  a  potential  calibration  equation.  Another  is 

obtained  by  subtraction  of  the  cold  reference  voltage  v  . 

o 

v  -  v  =  {rag'(0)]R  +  (r  -  r  )ag'(0).  (10) 

o  o 

Division  of  this  equation  by  the . amplification  factor  of  equation  (8) 
gives  yet  a  third  calibration  equation. 

(v  -  v  )/(v  -  v_)  -  [  x  /  ( r,  -  r  )  ]R  +  [(r  -  r  )/(r  -  r  )]  (11) 

oho  no  oho 

This  equation  expresses  the  output  voltage  as  a  fraction  of  the 

output  with  the  internal  reference  mirror  in  place,  all  voltages 

being  measured  above  the  cold  reference  level  v  . 

o 

Which  of  the  three  potential  calibration  equations  is  best? 

Equation  (11),  because  it  does  not  involve  the  response  factors  of 

the  detector  and  electronics  or,  more  precisely,  accounts  for  them 

through  v  and  v,  which  are  measured  before  each  scan.  Furthermore, 
oh 

the  behavior  of  the  slope  and  intercept  in  equation  (11)  can  be 

predicted;  this  will  provide  a  check  of  the  calibration's  success  or 

failure  and  perhaps  an  indication  of  what  went  wrong  in  the  event  of 

failure.  The  difference  (r.-^-r  )  will  increase  with  temperature 

h  o 

since  r,  includes  radiation  from  the  internal  reference  optics  and 
h 

the  radiometer's  interior.  However,  radiation  from  the  internal 

source  accounts  for  most  of  r,  ,  so  the  increase  of  ( r,  -  r  )  with 

h  ho 

temperature  should  be  slow.  This  will  result  in  a  slow  decrease  of 


the  slope  [t/(r^  -  r  ) )  with  temperature,  provided  the  loss  factor 

of  the  main  optics  depends  weakly  on  temperature.  The  difference 

(r  -  r  )  will  increase  with  temperature  more  rapidly  than  ( r,  -  r 
o  ho 

since  r  comes  entirely  from  the  radiometer's  interior  and  main 

optics.  The  intercept  [(r  -  r  )/(r,  -  r  )]  should,  therefore,  be 

o  n  o 

more  temperature  dependent  than  the  slope  and  should  increase  with 
temperature . 


NRL  FY85  BMAP  ACTIVITIES 


DATA  CALIBRATION 


Aug  84  Lab  Measurements  (Bedford) 

Much  learned 

MW  -  LARGE  OPTICAL  COUPLING 

LW  -  STABLE  SINCE  DEC  83  EMI  MEASUREMENTS 

Data  being  used  to  develop  cal  technique  &  software 

High  radiometer  temp  creates  stray  radiation  problems 

•  Questions 

Point  response  &  optical  coupling  vs  temp 
Radiometric  response  vs  temp 
Must  radiometer  temp  be  <  30  °C? 


Sept  84  Rooftop  Measurements  (Bedford) 

Accomplishments 

NRL  field  cal  protocol  looks  good 

540  calibrated  background  data  frames  being  disseminated 
Problems 

Intermittent  DAS  errors 

Slippage  of  internal  reference  outputs 


Jan  85  Lab  Measurement's  (Bedford) 

Point  &  extended  blackbody  measurements 
in  temp-controlled  environment 


Investigators:  Longmire,  Schaum,  Sweeney,  McHugh,  Lucke, 

Beaver 
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VUGRAPHS  AND  SUPPLEMENTAL  TEXT  FOR 
OPTICAL  SCIENCES  DIVISION  COLLOQUIUM,  17  JANUARY  1985 
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CALIBRATION 


RADIOMETER  DIAGRAM 


TYPES  OF  CALIBRATION  PERFORMED 


Point  Source 

(Point  Response/  Cross-channel  Coupling,  NED 

Bar  Source 

(Longwave  Resistive  Coupling  Correction) 


Near  Extended  Source 
(Radiometric  Response) 
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LABORATORY  CALIBRATION  PROCEDURE 


1.  Point  Source  Scans 

a.  Center  point-source  image  vertically  on  pixel. 

b.  Record  4  frames. 

C.  Do  A-B  FOR  ALL  CHANNELS  OF  DETECTOR  ARRAY. 

D.  Do  A-C  FOR  BOTH  ARRAYS. 

2.  Bar  Source  Scans 

a. '  Center  bar  target  vertically  on  detector  array. 

b.  Record  4  frames  in  alternating  mode  (2  frames 
from  each  detector) . 

c.  Do  a-b  with  bars  at  10.  20,  30,  40,  50  °C. 

.  3.  :  Near  Extended  Source  Scans 

a.  Remove  bar  apertures  from  source. 

b<  Place  source  close  to  telescope  objective  lens. 

c.  Record  4  frames. 

d.  Do  c  with  source  at  10,  20,  30,  40,  50  °C. 

4.  Do  1-4  with  radiometer  at  20,  25,  30,  35,  40  °C. 


1.  Point  Source  Scans 

.  a.  Center  image  from  portable  collimator 

VERTICALLY  ON  PIXEL. 

b.  Record  4  frames. 

c.  Do  a-b  for  channels  2.  8,  15  of  detector  array 

d.  Do  a-c  for  both  arrays. 

2.  No  BAR  SOURCE  SCANS. 

3.  Near  Extended  Source  Scans 

a.  Same  procedure  as  in  laboratory. 

b.  Record  scans  at  end  of  each  data  tape. 


1000  1500  2000  2500  3000 


SAMPLE  NUMBER 


Figure  4  —  Postamplifier  output  waveform  for  LWIR  channel  2 


Figure  5  —  Postamplifier  output  waveform  for  LWIR  channel  3 


SAMPLE  NUMBER 


Channel  3  of  corrected  data 


CALIBRATION 


ESULTS  AND  DISCUSSION 

Radiometric  response 


What  is  it? 

How  CAN  YOU  GET  IT? 


Answers  depend  on 

' INSTRUMENT 
SCENE 

FIELD  CONDITIONS 
DATA  USAGE 


l-uinu  wnvL  unu-iDr\n i  xuin 

TAPE  BP7 


(A-  A) /(  A-A)  =  0 


(dS « 7w/ 


TWO  APPROACHES  ARE  POSSIBLE. 

Derive  LW  response  function. 

Gorrect  LW  response  for  AC  coupling,  then  use  MW  response  equation. 


3  DAYS  9  HOURS  24  MINUTES  42,3  SECONDS 


CALIBRATION  STATISTICS  FOR  MID  WAVE  IR  DATA 


COEFFICIENTS 


CONFIDENCE 


3  CHANNEL 


IEL 

1 

• 

A  t  B 

3 

VARIANCE  3 

A  3 

B 

3 

1 

1 

l 

-0. 

195  :  2.  23 

a 

a 

3. 

591E-06  3 

2. 

319E-03  3 

2. 

927E-03 

a 

a 

2 

• 

1 

-0. 

201  1  2.  27 

"  I 

3. 

812E-06  3 

2. 

397E-03  3 

3. 

060E-03 

a 

a 

3 

1 

1 

-0.  232  !  2.  30 

• 

a 

5. 

194E-06  3 

2. 

B43E-03  3 

3. 

620E-03 

a 

a 

4 

• 

1 

-0. 

231  S  2.  32 

a 

a 

3. 

720E-06  3 

2. 

405E-03  3 

3. 

095E-03 

« 

a 

5 

a 

a 

-0. 

258  3  2.35 

3 

3. 

237E-06  3 

2. 

275E-03  3 

2. 

920E-03 

a 

a 

i) 

! 

-0. 

251  1  2.  36 

3 

7. 

109E-06  3 

3. 

359E-03  3 

4. 

36 IE— 03 

3 

7 

1 

a 

-0. 

268  S  2.  38 

a 

1. 

205E-06  3 

1. 

394E-03  3 

1. 

805E-03 

a 

• 

8 

• 

a 

-0. 

256  3  2.  38 

|  • 

4. 

916E-06  3 

2. 

800E-03  3 

3. 

651E-03 

a 

a 

9 

• 

-0. 

276  3  2.  41 

• 

a 

5. 

14 IE— 07  3 

9. 

14BE-04  3 

1. 

197E-03 

a 

a 

10 

• 

« 

-0. 

267  3  2.  43 

3 

5. 

771 E— 06  3 

3. 

051E-03  3 

4. 

037E-03 

a 

a 

11 

• 

• 

-0. 

273  3  2.  42 

a 

a 

5. 

943E-06  3 

3. 

105E-03  3 

4. 

081E-03 

a 

a 

12 

a 

• 

-O. 

250  3  2.  39 

a 

4. 

965E-06  5 

2. 

805E-03  3 

3. 

6B8E-03 

a 

a 

13 

a 

a 

-0. 

263  3  2.  40 

a 

a 

4. 

407E-06  3 

2. 

66 IE— 03  3 

3. 

487E-03 

a 

a 

14 

a 

a 

-O. 

244  3  2.  39 

a 

i 

4. 

624E-06  3 

2. 

699E-03  3 

3. 

557E-03 

a 

i 

15 

a 

a 

-O. 

251  3  2.  39 

a 

a 

4. 

152E-06  3 

2. 

567E-03  3 

3. 

367E-03 

a 

a 

16 

-i 

-O. 

227  3  2.  37 

t 

5. 

707E-06  3 

2. 

972E-03  3 

3. 

920E-03 

3 

PREDICTION  LIMITS 
FOR  MID  WAVE  IR  DATA 


3 

3  CHANNEL 

a  ' 
a 

a 

a 

0.  300 

a 

a 

0.  630 

U 

a 

a 

0. 

960  3 

1 

.  30 

> 

a 

a 

a 

3  1 

a 

a 

4.  271 E— 03 

a 

a 

4.  09 IE— 03 

T 

4. 

134E-03  3 

4. 

406E-03 

3 

3  2 

a 

• 

4.  399E-03 

a 

i 

4.  213E-03 

a 

a 

4. 

265E-03  ! 

4. 

559E-03 

a 

a 

3  3 

a 

i 

5.  145E-03 

a 

a 

4. 919E-03 

a 

» 

4. 

978E-03  3 

5. 

326E-03 

a 

a 

3  4 

a 

a 

4.  351E-03 

a 

« 

4.  161E-03 

a 

a 

4. 

218E-03  3 

4. 

524E-03 

a 

i 

3  5 

a 

» 

4.  066E-03 

a 

a 

3.  B83E-03 

a 

a 

3. 

934E-03  3 

4. 

223E-03 

a 

a 

3  6 

a 

a 

6.  019E-03 

a 

a 

5.  752E-03 

a 

a 

5. 

839E-03  ! 

6. 

282E-03 

a 

a 

3  7 

3 

2.  480E-03 

a 

a 

2.  368E-03 

3 

2. 

403E-03  3 

2. 

586E-03 

a 

a 

3  8 

» 

1 

5.  004E-03 

a 

a 

4.  783E-03 

a 

a 

4. 

859E-03  3 

5. 

234E-03 

a 

i 

3  9 

1 

a 

1 .  620E-03 

a 

a 

1.  547E-03 

a 

a 

1. 

572E-03  3 

1. 

697E-03 

a 

a 

3  10 

a 

a 

5.  42 IE— 03 

a 

a 

5.  180E-03 

a 

a 

5. 

275E-03  3 

5. 

707E-03 

a 

a 

3  .11 

i 

a 

5.  505E-03 

a 

a 

5.  258E-03 

a 

a 

5. 

349E-03  3 

5. 

779E-03 

a 

• 

S'  12 

a 

a 

5.  025E-03 

a 

a 

4.  805E-03 

a 

a 

4. 

898E-03  3 

5. 

274E-03 

a 

a 

3  13 

a 

i 

4.  739E-03 

a 

a 

4. 52BE-03 

a 

a 

4. 

604E-03  3 

4. 

968E-03 

a 

a 

3  14 

a 

i 

4.  847E-03 

a 

a 

4. 637E-03 

a 

a 

4. 

718E-03  3 

5. 

092E-03 

a 

a 

5  15 

a 

a 

4.  596E-03 

a 

a 

4.  395E-03 

a 

a 

4. 

469E-03  3 

4. 

B20E-03 

a 

a 

3  16 

a 

a 

5.  379E-03 

a 

a 

5.  151E-03 

a 

a 

5. 

241E-03  3 

5. 

652E-03 

a 

• 
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MIDWAVE  RADIOMETER  RESPONSE 

Monochromatic 

u  =  (v-v0)/(vH-v0) 

L  =  [(lh-l0)/t]u  -  [(l-l0)/t] 

Broadband 

<GTL>  =  <G(L.-Ln»U  -  <G(L-Ln)> 

HO  0 

=  Bu  +  A 

WHERE 

<Z>  =  /"Z’dX 

radiometer's  peak  normalized  spectral  response 
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Relation  between  radiometer  output  and 

BLACKBODY  RADIANCE  AT  RADIOMETER'S  ENTRANCE  APERTURE 


Requires  only 

BLACKBODY  STANDARD  SOURCE 

•..DETERMINATION  OF" SOURCE  RADIANCE  AT  ENTRANCE  APERTURE 
•  RADIOMETER'S  PEAK  NORMALIZED  SPECTRAL  RESPONSE 


Possible,  Calibration  equation  gives 

PEAK  NORMALIZED  BLACKBODY  RADIANCE 
WHICH  PRODUCES. SAME. RESPONSE  AS  PEAK  NORMALIZED  SCENE  RADIANCE 
WITH  BOTH  RADIANCES  MEASURED  AT  ENTRANCE  APERTURE. 


*F0R  DETAILS  SEE 

C.  L.  Wyatt,  Radiometric  Calibration  (Academic,  NY,  1978),  ch,  12. 


Relation  between  radiometer  output  and  scene  radiance 

AT  RADIOMETER'S  ENTRANCE  APERTURE 


Requires  standard  source  with 

SPECTRUM  AND  RADIANCE  RANGE 
OF  SCENE  RADIATION  AT  ENTRANCE  APERTURE 


Impossible.  Scene  radiation  at  entrance  aperture  varies  with 

SOLAR  ILLUMINATION 

meteorologic  conditions 

DISTANCE  TO  SCENE  OBJECTS 


BEST  CALIBRATION 


Relation  between  radiometer  output  and  scene  radiance 


Requires  standard  source 

WITH  SAME  SPECTRUM  AS  SCENE 
RADIANCE  RANGE  &  THAT  Of  SCENE 
PLACED  BESIDE  SCENE 


Impossible,  Scenes  contain  clouds 

WITH  VARIABLE  SPECTRA 
VARIABLE  RADIANCES 

AT  DIFFERENT  DISTANCES  FROM  RADIOMETER 


CALIBRATION  STATISTICS  FOR  LONG  WAVE  IR  DATA 


COEFFICIENTS 


CONFIDENCE 


CHANNEL 

:  a 

1 

B 

VARIANCE 

• 

1 

A 

B 

2 

-12.  9 

56.  2 

5.  109E-03 

t 

1 

0.  218 

0.  284 

3 

-12.  7 

56.  1 

5.  241E-03 

1 

1 

0.  220 

0.  287 

4 

-12.  7 

55.  9 

4.  030E-03 

• 

• 

0.  193 

0.  251 

5 

-12.  8 

56.  2 

3.  B48E-03 

1 

1 

0.  189 

0.  246 

6 

-12.  8 

55.  8 

3.  849E-03 

1 

• 

0.  189 

0.  245 

7 

-12.  5 

55.  6 

4. 435E-03 

1 

1 

0.  201 

0.  262 

8 

-12.  6 

55.  4 

4. 327E-03 

1 

> 

0.  199 

0.  257 

9 

-12.  8 

55.  5 

‘  5.  699E-03 

1 

a 

0.  230 

0.  296 

10 

-13.  0 

55.  4 

4.  379E-03 

• 

a 

0.  202 

0.  259 

11 

-12.  8 

55.  0 

4.  178E-03 

i 

a 

0.  197 

0.  251 

12 

-12.  8 

54.  5 

3.  277E-03 

a 

i 

0.  174 

0.  220 

13 

-13.  3 

54.  8 

5.  668E-03 

a 

i 

0.  232 

0.  291 

14 

-13.  2 

54.  0 

3.  852E-03 

a 

a 

0.  191 

0.  237 

15 

-14.  1 

54.  8 

5.  705E-03 

i 

a 

0.  237 

0.  292 

16 

-13.  7 

53.  9 

4.  439E-03 

i 

i 

0.  207 

0.  254 

PREDICTION  LIMITS 
FOR  LONG  WAVE  IR  DATA 


:  i  u 

i  CHANNEL  !  0.250  J  0.500  !  0.750  {  1.00 


2 

0.  213 

0.  174 

0.  1 58 

0.  173 

3 

0.  216 

0.  176 

0.  160 

0.  175 

4 

0.  189 

0.  154 

0.  140 

0.  153 

5 

0.  185 

0.  151 

0.  137 

0.  150 

6 

0.  185 

0.  151 

0.  137 

0.  149 

7 

0.  198 

0.  162 

0.  147 

0.  160 

8 

0.  196 

0.  160 

0.  145 

0.  158 

9 

0.  226 

0.  184 

0.  167 

0.  181 

10 

0.  199 

0.  162 

0.  146 

0.  158 

11 

0.  194 

0.  158 

0.  143 

0.  154 

•12 

0.  172 

•  0.  141 

0.  127 

0.  136 

13 

0.  228 

0.  186 

0.  167 

0.  178 

14 

0.  188 

0.  154 

0.  138 

0.  146 

15 

0.  231 

0.  188 

0.  168 

0.  177 

16 

0.  204 

0.  166 

0.  148 

0.  155 

CORRELATION: 


1.  00 


9  HOURS  29  MINUTES  3 1.8  . SECONDS 
LONG  WOVE  INFRA  RED 

Extended  source  T  =  50 . 0°C 


